NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



■ w 




' « 1 > 


■ ^ "V 


. r ( 




1- • 
it 

V 


tm 


. V, : ' ' ? 



\ )rc;\c 
l Jhiver sit\ 


% 


i 


College of Engineering 
Department of 
Mechanical Engineering 
and Mechanics 


PMUielphia, PA 19104 

(213)895-2352-53 


’ t 

» i 


(RAS A-CE- 176J73) CX1LES Oi- NITROGEN 
RfllSSluNj FBOfl 1h E CCfElSllCN CF 
NCNGUISPfchSE LltUiJ, FUEL SPRAYS Ph.D. 

Ihetixs (Diesel Uciv.) 173 p HC AOj/MF A01 

CSCL 13b G3/45 


Ntib- I 4770 


Unclus 

04930 


. ■ 

, £;■ 

- ' • tA ■ , ' %;r ‘ / 

• . • ' •; , A •• V • • •; "■ 





\ 




! « 




m,. 


OXIDES OF NITROGEN EMISSIONS FROM THE COMBUSTION OF 
MONODISPERSE LIQUID FUEL SPRAYS 


by 

Hamid Sarv 


Submitted in partial fulfillment of the requirements for 

the degree of 

Doctor of Philosophy 


Approved by 


7 ] '(? ' 


Dr. Nicholas P. Cemansjcy 
Principal Investigator 


Research Sponsored by the NASA Lewis Research Center 
through research grant NAG 3-1 


Department of Mechanical Engineering and Mechanics 
Drexel University, Philadelphia, PA 19104 


June 1985 


ABSTRACT 


A study of NO x formation in a one dimensional monodisperse spray 
combustion system, which allowed independent droplet size variation, was 
conducted. Temperature, NO and NO x concentrations were measured in the 
transition region, encompassing a 26-74 micron droplet, size range. Emission 
measurements of hydrocarbons, carbon monoxide, carbon dioxide and oxygen were 
also made. The equivalence ratio was varied between 0.8 and 1.2 for the fuels 
used in this study, including methanol, isopropanaol, n-heptane and n-octane. 
Pyridine and pyrrole were added to n-heptane as nitrogen-containing additives 
In order to simulate synthetic fuels. 

Results obtained from the postflame regions using the pure fuels indicate 
an optimum droplet size in the range of 43-58 microns for minimizing N0 X 
production. For the fuels examined, the maximum N0 X reductions relative to 
the small droplet size limit were about 10-20% for lean and 20-30% for 
stoichiometric and rich mixtures. This behavior is attributed to droplet 
interactions and the transition from diffusive to premixed type of burning. 
Preflame vaporization controls the gas phase stoichiometry which has a 
significant effect on the volume of the hot gases surrounding a fuel droplet, 
where N0 X is formed. On the other hand, the release of fuel vapor from the 
droplet to the reaction zone is influenced by droplet interactions which 
reduce the volume of the hot gases and N0 X production by bringing the flame 
closer to the droplet surface. 

For small droplets, gas phase combustion and droplet ignition are 
intensified by the high preflame extent of vaporization. Much of the N0 X is 
formed in the hot gases surrounding the burning fuel droplets, where premixed 
type of burning appears to dominate. Large droplets burn in a diffusive 
fashion, due to lower droplet vaporization before the flame front and the slow 
propagation of the fuel vapor emerging from the surface of the burning 
droplets. However, N0 X production for large droplets increases as they 
achieve elevated flame temperatures due to their higher Damkohler number. At 
some intermediate size, where minimum N0 X i s predicted, the flame temperature 


and the fuel vapor concentration In the gas phase fall in between the values 
observed for the small and large droplets. Under this situation, droplet 
interactions reduct the burning rate and lower N0 X formation by confining its 
production to a thin reaction zone. These interaction effects are also 
believed to be responsible for higher CO and O 2 , and lower CO 2 emissions. 

The occurrence of the minimum N0 X point at different droplet diameters 
for the various fuels appears to be governed by the fuel vaporization 
characteristics. This point was verified by air preheating and synthetic 
oxidizer experiments. 

Experiments with nitrogen-bearing additives revealed that under 
stoichiometric and rich conditions, reducing the droplet size increases the 
efficiency of fuel-N conversion to N0 X . This observation is associated with 
improved oxidation of the pyrolysis fragments of the additive by better oxygen 
penetration through the droplet flame zone. 

Preflame spray vaporization and N0 X formation from spray combustion are 
also studied analytically. Temperature and fuel, oxidizer and NO mass 
fractions around the burning droplets were determined from the numerical 
solution of the governing energy and species equations. A qualitative 
agreement between the predicted and experimental results was found. 
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ABSTRACT 


OXIDES OF NITROGEN EMISSIONS FROM THE COMBUSTION OF 
MONODISPERSE LIQUID FUEL SPRAYS 

Hamid Sarv 
and 

Nicholas P. Cernansky 

A study of NO^ formation in a one dimensional raonodisperse 
spray combustion system, which allowed independent droplet size 
variation, was conducted. Temperature, NO and NO concentrations 

X 

were measured in the transition region, encompassing a 26-74 
micron droplet size range. Emission measurements of hydrocarbons, 
carbon monoxide, carbon dioxide and oxygen were also made. The 
equivalence ratio was varied between 0.8 and 1.2 for the fuels 
used in this study, including methanol, isopropanol, n-heptane and 
n-octane. Pyridine and pyrrole were added to n-heptane as 
nitrogen-containing additives in order to simulate synthetic 
fuels. 

Results obtained from the postflame regions using the pure 
fuels indicate an optimum droplet size in the range of 43-58 
microns for minimizing N0 X production. For the fuels examined, 
the maximum N0 X reductions relative to the small droplet size 


limit were about 10-20% for lean and 20-30% for stoichiometric and 

rich mixtures. This behavior is attributed to droplet 

interactions and the transition from diffusive to premixed type of 

burning. Preflame vaporization controls the gas 1 phase 

stoichiometry which has a significant effect on the volume of the 

hot gases surrounding a fuel droplet, where N(^ is formed. On the 

other hand, the release of fuel vapor from the droplet to the 

reaction zone is influenced by droplet interactions which reduce 

the volume of the hot gases and NO production by bringing the 

x 

flame closer to the droplet surface. 

For oraall droplets, gas phase combustion and droplet Ignition 
are intensified by the high preflarae extent of vaporization. Much 
of the N0 x is formed in the hot gases surrounding the burning fuel 
droplets, where premixed type of burning appears to dominate. 

Large droplets burn in a diffusive fashion, due to lower droplet 
vaporization before the flame front and the slow propagation of 
the fuel vapor emerging from the surface of the burning droplets. 
However, NO production for large droplets increases as they 

X 

achieve elevated flame temperatures due to their higher Damkohler 
number. At some intermediate size, where minimum NO^ is 
predicted, the flame temperature and the fuel vapor concentration 
in the gas phase fall in between the values observed for the small 
and large droplets. Under this situation, droplet interactions 
reduce the burning rate and lower N0 x formation by confining its 
production to a thin reaction zone. These interaction effects are 
also believed to be responsible for higher CO and 0 2 , and lower C0 2 


emissions. 
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The occurrence of the minimum N0 X point at different droplet 
diameters for the various fuels' appears to be governed by the fuel 
vaporization characteristics. This point was verified by air 
preheating and synthetic oxidizer experiments. , 

Experiments with nitrogen-bearing additives revealed that 
un^r stoichiometric and rich conditions, reducing the droplet 
size increases the efficiency of fuel-N conversion to N0 X , This 
observation is associated with improved oxidation of the pyrolysis 
fragments of the additive by better oxygen penetration through the 
droplet flame zone. 

Preflame spray vaporization and N0 x formation from spray 
combustion were also studied analytically. Temperature and fuel, 
oxidizer and NO mass fractions around the burning droplets w;ere 
determined from the numerical solution of the governing energy and 
species equations. A qualitative agreement between the predicted 
and experimental results was found. 


CHAPTER 1 


Introduction 


Based on energy consumption figures published in 1982 
(Gustaferro, 1982), combustion of liquid hydrocarbon fuels 
constitutes 47% of the total energy usage in the United States. 
Despite the increasing contributions of other forms of energy such 
as nuclear, solar, and geothermal, our energy systems rely heavily 
upon the supply of fossil fuels. The forecast predicts that the 
combustion of liquid fuels will remain as the dominant source of 
power generation for years to come. Another intriguing fact is 
that about 42 percent of the liquid fuel is consumed in sprays. 

The ease in storage, handling, and atomization of liquid fuel has 
made spray combustion the driving force in large scale stationary 
gas turbines and boilers, as well as in air and some ground 
transportation vehicles. 

Because of the wide use of fuel spray combustion, the 
efficiency and pollutant formation from such applications have 
become of major concern to both government and private industry. 
Although there have been many analytical studies to characterize 
the combustion of fuel sprays, in most cases our knowledge is 
based on empirical relationships (obtained from experimental 
results), which are often scale dependent. As a consequence, 


there is still a lack of understanding of the most fundamental 
processes thal: occur during the combustion of fuel sprays* Most 
of the difficulty is associated with simultaneous occurrence of 
complex physical and chemical processes which make the formulation 
of a det-riptlve model of spray burning a strenuous task. 

Recently, a number of extensive reviews (Williams, 1973 ; Faeth, 
1979; Chigier, 1983; Sirigrmno , 1983) have discussed the state of 
the art in the area of liquid fuel spray combustion and have also 
suggested future directions for experimental and analytical 
developments. 

Oxides of nitrogen are among the major pollutants emitted from 

combustion sources. In spray combustion, NO (NO+NO ) emission 

X 2m 

levels have been found to be influenced by the fuel spray 
characteristics. While some researchers believe in a droplet size 
effect, others have attributed their observations to aerodynamic 
and mixing effects. Unfortunately, due to the operational 
inflexibilities in their experimental systems, the importance of 
the individual factors could not be distinguished. For instance, 
change of droplet size often required substantial variations in 
fuel and air flow rates which may have influenced the combustion 

t 

processes. Also, the related analytical studies have deviated 
from modeling the physical system by limiting their attention to 
the combustion of single droplets. This study was undertaken to 
resolve some of the contradictions by isolating the important 
parameters that could lead to lower NO formation. 

A well-characterized spray combustion facility was used to 
simulate some of the basic processes that occur in practical 


combustors such as diesel engines or gas turbines. The 
experimental variables include droplet size* stoichiometry* fuel 
type, oxidizer composition* etc. Oxides of nitrogen and 
temperature constitute the major portion or* the experimental 
measurements made in this study. 

Experimental results that are presented for a range of 
operating conditions indicate an optimum droplet size for 
minimizing NO production. Droplet interactions and transition 
from diffusive to premixed type of burning are considered to be 
responsible for the observed behavior. This point has also been 
examined through additional measurements and analytical 
calculations. Numerical predictions are compared with the 
experimental results and possible refinements to the model are 
discussed. 

Chapter 2 includes a review of some of the previous work in 
the area of spray combustion and pollutant formation. The 
experimental facility and measurements are reported in Chapters 3 
and 4, respectively. The analytical calculations and predictions 
are presented in Chapter 5. Finally, the conclusions of the study 
and recommendation for future work are given in Chapter 6. 


CHAPTER 2 


Background and Literature Review 

As noted, a number of extensive reviews dealing with fuel 
spray combustion have appeared in the literature (Williams, 1973; 
Faeth, 1979; Chigier, 1983; Sirignano, 1983). No attempt is made 
to duplicate or summarize these landmark efforts here. Instead, 
only the previous work particularly germane to the present study 
is reviewed in this chapter. The three major areas covered are 
liquid fuel spray combustion, transition region phenomena and N0 X 
formation. Some of the specific topics such as droplet 
interactions and NO^ formation mechanisms for fuels with 
chemically bound nitrogen are examined more closely along with the 
relevant findings of this study, in later chapters. 

2.1 Liquid Fuel Spray Combustion 

The most fundamental aspect of spray combustion involves the 
burning of single fuel droplets. Experiments with droplets 
suspended from filaments (Godsave, 1953), porous spheres fed by 
liquid fuel (Spalding, 1953), freely falling droplets (Kumagai and 
Isoda, 1956), and dilute sprays (Bolt and Boyle, 1956) have shown 
that the surface regression rate of a fuel droplet during 
evaporation or combustion varies linearly with time and is 


2 

governed by the ”D law”. For Isolated droplets, theoretical 
calculations relate this rate to the droplet size, fuel/ambient 
gas composition, local temperature, pressure, and the relative 
velocity between the fuel droplet and the gas phase. Excellent 
agreement has been found between the physical model and the 
experimental results. 

Development of a chemical model for predicting the 
concentration of the combustion products requires consideration of 
heat, mass, and momentum transfer in addition to a knowledge of 
the major chemical reactions. Certain assumptions, including 
single fuel droplet, spherical symmetry, equal thermal and 
molecular diffusivitles, and a one step fuel oxidation reaction, 
are commonly made to simplify the analysis. In spite of the 
simplifications, the problem remains complex and numerical 
treatment becomes necessary. One of the first analyses of this 
type was done by Lorell et al. (1956) who determined the flame 
structure around a single ethanol droplet under steady state 
conditions. In their analysis, the flame temperature and its 
thickness were found to be sensitive to the chemical reaction 
rates. When the reactions become infinitely fast, the reaction 

t 

zone reduces to a flame shell of infinitesimal thickness. In this 
situation, temperature and concentration profiles as well as flame 
location can be determined explicitly (Williams, 1965). Unsteady 
combustion of a single fuel droplet placed in a hot and reactive 
environment was numerically investigated by Botros et al . (1980). 
Their study showed that dif fusional or premixed type of burning : 
depends on the fuel concentration in the droplet environment. In 


2 

contrast with the D law predictions, they demonstrated an 
unsteady flame front relative to the receding droplet radius. 

New complications arise when one attempts to characterize the 
behavior of fuel sprays. For instance, droplet size distribution 
and air flow patterns require special consider ^.ons due to their 
effects on spray combustion (Beer and Chigier, 1983). Spray 
vaporization results in fuel enrichment and cooling of the gas 
phase (Law, 1977; Rao and Lefebvre, 1981) which can have 
significant effects on ignition and burning of the droplets. 

An individual droplet in a spray has to compete with its 
neighbors for the available oxygen. This competition becomes 
stronger for smaller droplet spacings and leads to an increase in 
the droplet lifetime compared to that of an isolated droplet with 
the same diameter. 

Kanevsky (1956) showed that for a stationary array of nine 
droplets, the burning rate for the center droplet reduced when the 
droplet flame boundaries merged at decreased droplet spacings. 
Similar results have also been obtained by others (Nuruzzaraan et 
Si* , 1971; Miyasaka an^ Law; 1981; Sangiovanni and Labowsky, 1982; 
Xiong et^ al . , 1985) using linear arrays of monodisperse fuel 
droplets. 

The theoretical work of Labowsky ( 1976-1980) , Umemura et al . 
(1981a, 1981b) , and Marberry et al. (1984), assuming a quasi-steady 
mass flux transport model, show that the evaporation rates .of 
interacting droplets decrease when the droplet spacing is reduced. 
Their calculations indicate that the spray evaporation rate can be 
determined by applying a correction factor to the single droplet 


model which Is only dependent on the array geometry. 

Chiu gi al. (1977,1983) have Identified four classes of sprays 
that exhibit external sheath, external group, internal group, and 
single droplet modes of combustion. These regimes are separated 
by a group combustion number which is defined as the ratio of the 
heat transfer in the gas phase to the heat transfer between the 
gas and the liquid phase. In a recent study, Bellan and Cuff el 
(19 c ;3) presented a spray evaporation model in which the 
interacting droplets were separated by spheres of influence, and 
each droplet in the spray evaporated inside a finite surroundings. 

Their predicted results show that the evaporation of interacting 

2 

droplets cannot be described by the D law over the entire droplet 
lifetime. 

2.2 Transition Region Effects 

Another point of interest in spray combustion is the behavior 
of spray systems in the "transition region". In this region, 
which covers a range of 15-80 pm droplet diameter, the local 
conditions around the droplets are influenced by their size and 
spacings. Fcr these sprays, the dominance of heterogeneous or 
homogeneous type of burning depends on the amount of 
prevaporization before the flame and * '.a strength of droplet 
interactions in the flame. Such factors have been determined to 
be responsible for a number of interesting and important 
observations including increased burning velocities 
(Polymeropoulos and Das, 1975; Hayashi et al. , 1976; 
Polymeropoulos, 1984), broadened flammability limits (Burgoyne and 


Cohen, 1954; Hayashi et al., 1981), lower Ignition energies (Chan 

and Polymeropoulos, 1981), and reduced NO emissions (Nizami et 

x 

al, , 1978-82). These findings are of significant value for the 
design and development of high efficiency and low emission spray 
combustors. 

2.3 NO Formation in Spray Combustion Systems 

X 

Prediction of NO emissions from spray combustion sources by 
isolating the contributing factors has been the subject of many 
investigations. Nitric oxide formation in the vicinity of a 
diffusion flame around a single droplet has been analytically 
treated by Kesten (1972) and Bracco (1973). The former study 
dealt with transient NO formation, assuming quasi-steady droplet 
combustion and an infinitely fast fuel oxidation reaction, while 
in the latter case, steady state droplet combustion at a finite 
reaction rate was considered. Both studies showed that the NO 
production per mass of fuel (Emission Index) increases with 
increasing droplet diameter. Kesten also demonstrated that under 
certain constraints (residence time equal to the droplet lifetime 
and near stoichiometric values), the NO emission index for 
homogeneous gas phase combustion can be more significant than that 
from the droplet combustion. 

Exit plane measurements from a variety of experimental and 
practical combustion systems have also indicated the effects of 
droplet size and spray characteristics on NO emissions. 

X 

Increased NO emissions with increasing Sauter Mean Diameter (SMD) 
have been measured from a marine gas turbine engine (Opdyke, 


1983). Mellor and co-workers (1976, 1980) showed that the N0^ 
emissions from a gas turbine combustor can be correlated to the 
characteristic times for droplet evaporation, mixing, and chemical 
reactions. In a different study using a research gas turbine 
combustor, Nicolls jet al . (1980) reported a decrease in the NO 
levels when the droplet size was increased from 10 to 57 microns 
at a constant residence time. 

Murphy and Borman (1979) evaluated the contrioution of droplet 
burning to the total N0 x formed in their laminar mist flame. 

Their result indicated higher NO^ emissions for larger droplets. 
With a monodisperse fuel spray burning in an aerodynamically 
stabilized flame inside a combustion tunnel, Nizami and Cernansky 
(1978-1982) showed the existence of an optimum droplet diameter 
(around 50 pm) for minimum NO production. 

Combustion of fuels which contain chemically bound nitrogen 
contributes to additional NO emissions over what is formed from 
atmospheric nitrogen. Unfortunately, very little is known about 
the effect of droplet diameter on fuel-N conversion to N0 x . 

Perhaps the only qualitative study in this area was done by 
Appleton and Heywood (1973) who showed a strong influence of the 
degree of atomization and fuel/air mixing on the conversion 
efficiency. In a similar study, Foster and Keck (1980) reported 
significant levels of hydrogen cyanide in fuel-rich turbulent 
diffusion flames as far as five burner diameters downstream of a 
fuel atomizer nozzle. Increasing atomization pressure led to 
intense mixing and resulted in major reductions in HCN and a 
moderate rise in NO concentrations. Crowhurst and Simmons (1985a) 


utilized a hollow sintered cylinder burner for a structural 
examination of methane-acetonitrile-argon diffusion flames burning 
in a stream of oxygen and argon. They concluded that nitric oxide 
forms primarily on the oxygen side of the main reaction zone, with 
some of the NO diffusing back into the pyrolysis region of the 
flame and reacting with the nitrogenous radicals to form molecular 
nitrogen. 

Most mechanistic studies have concentrated on prevaporized and 
premixed fuel/air mixtures doped with a nitrogen containing 
compound. Overall, the N0 X formation reactions from fuel-N 
involve the oxidation of the pyrolysis fragments of the nitrogen 
additive. Such mechanisms are affected by the flame stoichiometry 
and have been discussed by Feniraore (1976a, 1976b), Haynes 
(1977a), Crowhurst and Simmons (1983, 1985a, 1985b), Dasch and 
Blint (1984), and others. Previous work of Fenimore (1976a), 
Foster and Keck (1980), Sapre and Quader (1983) and Chen and Malte 
(1984), have shown little or no effect of the type of doping 
compound on N0 X yield. 


CHAPTER 3 


Experimental Set-Up and Procedures 

3.1 Introduction 

In this chapter the basic combustion facility and the N0 X and 
temperature measurement systems are described. The operating 
conditions and the actual measurements are discussed in Chapter 4. 
Other instrumentation including the exhaust analyzers and the fuel 
vaporizing system are specified in Appendix A. 

3.2 Spray Combustion Facility 

3.2.1 Monodisperse Spray Generation System 
A Berglund-Liu Vibrating Orifice Monodisperse Aerosol 
Generator, Model 3050 (Berglund and Liu, 1973), is used to supply 
a monosized spray of droplets to the combustion facility. This 
generator has been employed successfully in other monodisperse 
spray combustion studies (Nizami et al . , 1978-1982; Koshland and 
Bowman, 1983, 1985). 

Monodisperse droplet generation is achieved by applying 
periodic disturbances to a piezoelectric ceramic which exerts 
mechanical vibrations to a liquid jet emerging from an orifice 
plate. The oscillation frequency is applied from a function 
generator (Hewlett Packard Model 3310A). The perturbed jet then 


breaks up Into discrete droplets with a standard deviation of 

approximately 1% of the mean diameter. Since one droplet is 

generated per cycle cf disturbance, the droplet diameter can, be 

calculated iron the liquid feed rate (Q) and the, frequency of 

1/ 3 

disturbance (f) and is given by : D*(6Q/nf) . Thus, droplet size 
can be varied by changing the frequency or the fuel flow rate. 

Although the actual orifice size is not used in the 
calculation of droplet diameter, an orifice plate is only capable 
of providing monosized droplets within a relatively narrow range. 
However, the monodispersion range can be extended easily by using 
a set of orifice plates with various orifice diameters. 

3.2.2 Fuel and Air Flow Systems 

A schematic of the experimental facility is shown in Figure 
3.1. A high infusion syringe pump (Harvard Apparatus Model 901) 
supplies liquid fuel to the aerosol generator. Air is supplied at 
a regulated pressure of 35 psig from the house air compressor. 

The stream of uniform droplets emerging from the orifice plate 
is subjected to a turbulent jet of dispersion air to prevent 
coagulation, A secondary flow of dilution air is also added to 
achieve the desired stoichiometry. For minimum fluctuations in 
the flow rates, the air flows are adjusted by electronic flow 
controllers (Tylan Model FC-260), and monitored by Mass Flow 
meters (Hastings Model ALL-5K). 

A tapered baffle plate with 3 cm I.D. is positioned at 2 cm 
above the droplet generator to improve the mixing of the dilution 
air with the dispersed aerosol. The air/fuel aerosol then passes 








through a transition piece with smooth curvatures designed to 
minimize any recirculation. Above the transition piece are flow 
reducing sections which direct the laminarlzed aerosol to a flame 
holder where the monodisperse fuel spray is burned inside a screen 
stabilized one dimensional flame. The flame is surrounded by a 
long Pyrex tube to avoid outside air entrainment. 

3.2.3 Flame Holders 

In order to maintain one dimensional and stable flames for all 
equivalence ratios, a series of 40x40 mesh stainless steel screen 
flame holders (0.25 mm wire diameter) with matching flow reducing 
sections were designed and fabricated. The wire meshes block 50% 
of the area. In the absence of a flame, the fraction of droplets 
that survive the screen flame holder was determined using a dual 
beam He-Ne laser operated in the back scatter mode. By traversing 
the focal volume across the burner, a count rate of the droplets 
was determined both with and without the screen in place. The 
measurements indicate that the percentage of droplets which 
survive impaction on the screen are 37 and 24% for the 37 and 70 
pm droplets, respectively. In the presence of a flame, those 
droplets which do not make it through the screen and agglomerate 
on the screen wires will flash evaporate and add to the vapor 
phase in the fuel aerosol entering the flame. The flame holders 
have exit diameters of 9.5, 12.7, 14.3, and 17 mm and are used for 
equivalence ratios between 0.8 and 1.2. Both the flame holder and 
the flow reducing section are water cooled to prevent uncontrolled 
droplet vaporization. 


3.2.4 Other Design Considerations 

The design of the spray burner facility allows for adjustment 
in the burner height from 16-23 era above the orifice plate. A 
burner height of 17.8 cm was selected experimentally based on the 
aerodynamics of the dispersed aerosol. Below this optimum burner 
height, the fuel spray jet is not completely spread and diffused* 
resulting in a nonuniform velocity and fuel/air concentration at 
the screen* as indicated by the shape and luminosity of the flame. 
Above this optimum height, a significant number of droplets begin 
to settle out and accumulate on the droplet generator, leaving 
them unburnt. Thus, the burner height is fixed by these practical 
considerations. 

3.3 Sampling and Analytical Instrumentation 

3.3.1 N0 X Measurement System 

Nitric oxide and total N0 X are collected using an uncooled 
aerodynamic stainless steel sampling probe with an orifice 
diameter of 0.3 mm as shown in Figure 3.2. For comparative 
studies of the effect of sampling probes on N0 X measurements, a 
water cooled probe having an orifice size of 0.25 mm (Figure 3.3) 
is used. The probe is cooled to 50 *0 during all measurements. 

The gases are sampled at a line pressure of 80 torr to prevent 
water condensation in the sampling line. Initially, the probes 
were designed to achieve aerodynamic quenching. However, Colket 
et ed. (1982) have argued that a single probe may not be capable 
of aerodynamic quenching for all ambient conditions. But, as 





pointed out In section 4.7.1, the use of cooled or uncooled 
sampling probes did not have a noticeable effect on the measured 
N0 X values and only resulted in differences between the NOj 
levels. 

The concentrations of NO and N0 X are measured by a 
chemiluminescent NO analyzer (TECO Model 12A), equipped with a 
Molybdenum NO^ to NO converter. This system had been modified for 
low pressure sampling (Nizami, 1978). A high vacuum sampling pump 
and a variable sample capillary are used to maintain the analyzer 
reaction chamber pressure at 8 torr. For each steady state burner 
operating condition, continuous N0 X sampling and measurements are 
carried out for a period of two minutes and the measurements are 
plotted by a strip chart recorder for accurate N0 X concentration 
readings. Using this method, all N0 X readings were accurate to 
within 1 ppm, with a day to day reproducibility of 5 percent. 

3.3.2 Temperature Measurement System 

Temperature measurements are made with a Pt/Pt-13% Rh fine 
wire thermocouple (0.075 mm diameter) with a bead diameter of 
approximately 0.11 mm. The wires are inserted into a two-hole 
ceramic support which is attached to the side of the uncooled 
sampling probe (see Figure 3.2). This arrangement allows 
simultaneous N0 X and temperature measurements. The thermocouple 
is connected to a digital display Trendicator (Doric Type 400 A) 
as well as to a strip chart recorder for accurate measurements of 
the mean temperature at a particular location in the burner 
system. Positioning of the combination gas sampling/thermocouple 


probe Is achieved by a 3-D traversing mechanism. 

In order to determine the actual temperature at a sampling 
position, the measured value must be corrected for conduction and 
radiation losses. Such corrections are laborious and contain some 
level of uncertainty as discussed in Anpendix D, For this reason, 
no corrections were applied and the temperature results reported 
In Chapter 4 represent the raw measurements only. Uncorrected 
temperature readings are accurate to within 20 *C, with a day to 
day reproducibility of 5 percent. 


CHAPTER 4 


Experimental Results and Discussion 

4.1 Introduction 

Selected results on the effects of fuel type, stoichiometry, 
air preheating, oxidizer composition and fuel nitrogen additives 
on temperature and NO^ formation are presented and discussed in 
this chapter. Additional experimental results are given in 
Appendix B, 

In some cases, slight modifications to the experimental 
facility were necessary to attain the desired operating 
conditions. Those changes are described appropriately. 

The burner was operated under atmospheric pressure and in the 
monodisperse droplet size range at a constant fuel feed rate of 
0.382 cc/min. To achieve the desired stoichiometry, the dilution 
air flow rate was varied, keeping the dispersion air flow rate at, 
1300 cc/min, 

4.2 System Verification 

4.2.1 Monodispersion Range 

Since a single orifice plate can only operate in a narrow 
monodisperse range of droplets, a set of 20, 25, and 35 pm orifice 
plates were selected to cover a combined range of 36 to 74 pm 


droplet diameter. Figure 4.1 shows the raonodisperse operating 
range of the aerosol generator determined by the manufacturer's 
suggested "jet deflection" method (Thermo System Inc., 1977) and 
further confirmed by a slide impaction technique recommended by 
Berglund and Liu (1973). All droplets in a collected sample must 
have the same size in order to be considered monodisperse. 

Typical slide impaction results of raonodisperse droplets are shown 
in Figure 4.2. 

4.2.2 Radial Temperature and NO^ Mapping 

One dimensionality and radial symmetry of the flame for three 
different droplet sizes were examined by temperature and N0 X 
traverses at a height of 4 mm across the 17 mm flame holder. 

In general, the flames appeared as pale blue and, for small 
droplets > it looked as if tiny streaks were burning in a 1 ram 
thick envelope flame with a luminous region exfcented to about 3 
ram. Large droplets seemed to have occasional merging individual 
flaraelets and ..ere luminous up to 6 mm above the flame holder. 

The luminous region was relatively longer for lower volatility 
fuels. 

The radial temperature and NO* profiles shown in Figures 4.3 

* 

and 4.4 indicate some asymmetry for the 69.1 pm droplets, but for 
most part they are flat, thus confirming the one dimensionality of 
the burner system. 


* A 15 pm orifice plate was used to produce the 26.7 pm droplets. 





37.5 micron droplets 70.0 micron droplets 



figure 4.2. Microscopic photographs of droplet impactions 
glass slides. 
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4.3 Fuel Type and Stoichiometry 

Post flame temperature, NO, and N0 X were measured in the 36 to 
70 pm droplet size range. Several hydrocarbon fuels were used 
including isopropanol, methanol, n-heptane, and n-octane. All 
fuels, except for n-octane were HPLC grade and had a purity of 
greater than 99 percent (n-octane had a purity of 95%). A fixed 
sampling location on the burner centerline at 5 cm above the 
screen flame holder was selected to avoid regions of high 
temperature and concentration gradients found in axial probing. 
These results are shown in Figures 4. 5-4. 8. Where available, 
prevaporized/premixed limit results are also shown as data points 
at zero droplet diameter. For these premixed experiments, liquid 
fuel was vaporized in an evaporator and the fuel vapor was mixed 
with the dispersion air before entering the combustion facility. 

In all caseg, a general trend was observed in which both NO and N0 X 
decreased with reducing droplet diameter, reaching minima around 
43-58 pm, and then increasing to fairly constant levels with any 
further decrease in the droplet diameter. The equivalence ratio 
and fuel type appeared to only affect the absolute levels of N0 X 
measured, without appreciably affecting the observed droplet 
behavior. However, N0 X reductions were larger for rich mixtures 

with increased number densities and droplet interactions. For 

0 

these tests at an air temperature of 24 C, the N0 X minima occurred 
at 43, 49* 54, and 58 pm initial droplet diameters for n-octane, 
n-heptane, isopropanol, and methanol, respectively. 
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Figure 4.6. Droplet size effect on post flame temperature, NO and 
NO*, for methanol at various equivalence ratios, 
o - NO; • - NO^; A - Temperature. 
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Figure 4.7. Droplet size effect on post flame temperature, NO and 
NO*, for n-heptane at various equivalence ratios, 
o - NO; • - N0 X ; A- Temperature. 
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4.8. Droplet size effect on post flame temperature, NO and 
N0 X , for n-octane at various equivalence ratios, 
o - NO: • - NO*: A- Temperature. 


Similar trends in the experimental data were also £ound by 
Nizami et al . (1978-82) who hypothesized the following mechanism 
to explain the observed behavior. Decreasing the droplet diameter 
at a constant fuel feed rate Increases the droplet number density 
and, hence, droplet interactions, leading to depletion of the 
local oxygen concentration around the burning droplets and a 
suppression of the local flame temperature and NO* formation, 

This analysis is also consistent with the results of Nuruzzaman et 
al. (1971), Kesten (1972), and Bracco (1973). Continued reduction 
in droplet size further Increases number density, but enhanced 
vaporization brings about increased droplet spacings, mitigates 
droplet interactions, increases burning rates and flame 
temperature, and eventually reverses the decreasing N0 X trend as 
the situation approaches the premixed case. 

The variations in N0 X minima appear to be governed by the 
differences in the physical properties of the individual fuels. 
Larger initial droplet diameters are needed for high volatility 
fuels, such as methanol, to overcome rapid prevaporization and to 
give some level of droplet interactions in the flame. By the same 
argument, lower volatility fuels like n-octane, can achieve the 
optimum conditions at smaller droplet diameters. 

In the next two sections, additional experimental results on 
the effect of extent of prevaporization and droplet interactions 
are presented. These phenomena have been incorporated into an 
analytical model that is described in Chapter 5. 


4.4 Air Preheating and Prevaporization 

In the previous section the shifts In the N0 X minimum points 
for the test fuels were found to be consistent with the 
differences in their evaporation and combustion characteristics. 

a 

In this section, the importance of enhanced prevaporization in 
addition to droplet size effects on N0 X formation were 
investigated by air preheating experiments. 

Air preheating was accomplished by wrapping heating tapes 
around the upper and lower housing assemblies (Figure 3.1) and 
controlling the temperature by a variac. For a constant setting, 
there were only small variations in temperature (less than one 
degree centigrade) throughout the combustor. The water flow rate 
in the cooling jacket was also adjusted to keep the exit water 
temperature from the jacket equal to that of the preheated air. 

For this study, measurements were made at a height of 10 cm above 
the screen flame holder (the adjustment in sampling location was 
made to reduce potential catalytic NOg-* NO conversion reactions 
in the uncooled sampling probe). 

The effect of air preheating and initial droplet diameter on 
temperature, NO, and N0 X for the various fuels are shown in 
Figures 4.9-4.12. Air preheating results in enhanced vaporization 
of droplets, shifting the minimum N0 X point towards larger initial 
droplet sizes and mitigating the reduction in N0 X levels at the NO 
minimum point. These effects are believed to result from a 
balancing of two phenomena: (1) a requirement for larger initial 
droplet sizes to maintain some optimum extent of prevaporization, 
since air preheating increases prevaporization and local fuel 
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Figure 4.9, Air preheating and droplet size effects on post flame 
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vapor concentrations as well as droplet spacings; and (2) fewer 
droplets and reduced droplet interactions as the size is 
increased. This tradeoff is particularly pronounced in the case 
of methanol (Figure 4.9), where the N0 X reduction phenomena almost 
disappears with preheating. 

In order to ad&qn&tely interpret these results, it is 
important to realize that the effectiveness of air preheating on 
fuel vaporization relies upon the droplet travel time as well as 
the specific fuel properties. At a constant equivalence ratio and 
fuel feed rate, the air flow requirement increases with increasing 
fuel molecular weight, which consequently reduces the droplet 
travel time. For instance, N0 K minima for methanol, having the 
highest volatility and droplet travel time among the test fuels, 
quickly moved towards larger initial droplet diameters and the N0 X 
levels began to reach a constant value at elevated air 
temperatures. These effects were somewhat retarded for 
isopropanol due to its lower volatility and shorter droplet 
residence time. 

For n-heptane, there is a relatively shorter droplet residence 
time. Thus, preheat temperature increases primarily affected the 
small droplet sizes and their associated NO* values, with some 
smoothing in the large droplet region. N-octane, with the lowest 
volatility and shortest droplet travel time, seemed relatively 
unaffected by the increases in temperature. Air preheating 
temperatures above 44 °C are required to cause a significant shift 
in N0 X minima for n-octane. These factors are discussed in more 
detail as part of the analytical studies in Chapter 5. 


4.5 Synthetic Oxidizers 


Results of the parametric fuel studies and air preheating 
experiments indicated the importance of parameters such as fuel 
properties and temperature on droplet diameter and interactions. 

In spray combustion, changes in flame temperature affect the 
droplet lifetime, interactions with neighboring droplets, and the 
formation of NO*. These effects were studied by a series of 
experiments in which air was partly replaced by a mixture of 
oxygen and argon, keeping all other operating conditions the same. 

Air, oxygen, and argon were individually metered, then mixed 
in a cross, and later divided into the usual dispersion/dilution 
proportions before entering the combustion facility. In all cases 
the volumetric ratio of oxygen to the inert gas (nitrogen plus 
argon) was kept at 2X/79 . 

The experimental results shown in Figure 4,13 demonstrate the 
effect of droplet size and fraction of synthetic oxidizer on NO, 
N0 X , and temperature at a sampling height of 10 cm above the flame 
holder. Since argon is a monatomic gas with a lower heat capacity 
than nitrogen, its substitution for nitrogen increases the flame 
temperature and flame speed. For this reason, different size 
flame holders had to be used (as indicated on the Figure) to 
stabilize the flame. 

Compared to air, the N0 X minimum point moves towards larger 
initial droplet diameters as the fraction of the synthetic 
oxidizer increases in the gas phase. Lesser N0 X reductions and 
higher N0 X levels (relative to air) can also be observed from the 







All of the above observations are due to changes in flame 
temperature. An increase in temperature reduces the droplet 
lifetime and diminishes the droplet interactions in flame. As a 
consequence, the minimum N0 X point occurs with larger droplets 
which burn for a longer time. However, shifts towards larger 
droplet diameters are accompanied by a lower droplet number 
density and reduced interactions. Therefore, at very high 
temperatures, the effect of droplet interactions on N0 X production 
is minimal. 

4.6 Fuel-Nitrogen Additives 

A series of experiments with simulated synthetic fuels were 
conducted in order to investigate the effect of droplet size on 
the conversion of fuel-nitrogen to N0 X . In general, synthetic 
fuels contain traces of nitrogen-containing organic compounds 
including pyridines, nitriles, quinolines, pyrroles, indoles and 
carbazoles (Chan jat al. , 1984) . These compounds undergo thermal 
decomposition in the reaction zone, producing radicals such as 
HCN, N, CN, and NH which can either oxidize to NO* or recombine to 
form molecular nitrogen. 

N-heptane was doped with pyridine (C^H^N) and pyrrole (C^H^NH) 
to give a 0. 1-1.0 % nitrogen content by weight in the fuel 
(typically found in synthetic fuels from shale oil). 

The experimental data collected at a sampling position of 10 
cm above the burner (Figures 4.14 and 4.15), indicate that the 
optimum droplet diameter for minimizing NO^ (which was observed 
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previously with pure n-heptane) disappears with Increases in the 
fuel-N concentration. Also, it was noticed that more fuel-N is 
converted to NO^ at small droplet diameters. 

Appleton and Heywood (1973) suggested that the degree of 
air/fuel mixing and the availability of oxygen atoms control the 
oxidation of the nitrogenous species. In their investigation, 
maximum conversion efficiencies were found for fuel-lean and well 
atomized cases. The present results show that under 
stoichiometric and rich conditions, reducing the droplet size 
improves the oxidation of the pyrolysis fragments of the additives 
to NO* by better oxygen penetration through the droplet flame 
zones. The abundance of oxygen for lean mixtures made this 
process less dependent on the size of droplets as shown in Figures 
4.14 and 4.15. Changing the dopant from pyridine to pyrrole with 
a different molecular structure did not change the results. The 
fact that the type of nitrogen-bearing additive does not change 
the N0 X yield has also been addressed in a number of other studies 
(Fenimore, 1976a; Foster and Keck, 1980; Sapre and Quader, 1983). 

It should be pointed out that the presented results include 
the thermal N0 X contribution. However, for fuel-N, the dependence 
of nitric oxide yield on flame temperature is not as strong as it 
is for the oxidation of molecular nitrogen (De Soete, 1974; 
Fenimore, 1976a). Therefore, temperature changes resulting from 
droplet interactions become unimportant as the concentration of 
the nitrogen-containing compounds in the fuel increases. 

The conversion efficiency of the fuel nitrogen to NOy 
(determined from a nitrogen balance and subtracting the thermal NO 


contribution) for the 37 pm droplets are given in Table 4.1. The 
fraction of fuel-N that was not recovered as NO^ presumably 
appeared mostly as Ng in the products. Foster and Keck (1980) 
measured high HCN concentrations in fuel-rich turbulent diffusion 
flames as far as five burner diameters downstream of the fuel 
nozzle, constituting a large portion of the nitrogenous species 

(hcn+no+nh 3 ) . 

In the following section, inflame and post flame measurements 
of NO, N0 X , and temperature are presented. 

4.7 N0 X Structure in Flame 

4.7.1 Pure (Nitrogen Free) Fuels 

Axial temperature, NO and NO* concentrations were measured for 
various fuels and droplet sizes by traversing the gas sampling/ 
thermocouple probe along the burner centerline inside a 27 cm long 
Pyrex tube. Figure 4.16-4.19 represent the measurements for 
stoichiometric combustion of methanol, isopropanol, n-heptane, and 
n-octane aerosols, respectively. In general, temperature, NO and 
N0 X concentrations increase rapidly in the flame zone, with NO and 
NO* reaching their maximum values slightly above the temperature 
peak and then, decreasing gradually to relatively constant levels. 

The effect of sampling probe on the measured NO and NO* values 
was also studied using the water cooled probe. The probe was 
cooled to 50 *C during all measurements. Figures 4.20-4.22 show 
the concentration profiles of NO and NO* sampled by both cooled 
and uncooled probes. The results indicate minor differences 
between the NO* values for each condition while a higher NO^ yield 


* 

Table 4.1. Conversion efficiency of fuel-N to NO* from 
the combustion of 37 pm n-heptane droplets 
doped with pyridine. 


% N By Weight of Fuel 

<f)*>0.85 

<t>-1.0 

<j)-1.2 

0.1 

73 

74 

69 


(70) 

(69) 

(55) 

0.5 

70 

80 

59 


(67) 

(70) 

(48) 

1.0 

70 

76 

59 


(67) 

(66) 

(48) 


* Calculations are based on the maximum NO* concentration measured 
along the burner centerline. Values inside the parentheses are 
calculated based on the steady state NO^ concentration at 10 cm 
above the burner. 
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is observed when gases were sampled with the cooled probe «, This 
occurs as a result of the oxidation of NO by HQ^ which is formed 
under rapid quenching in the sampling probe (Cernansky, 1977; 

Hori, 1980; Duterque et al., 1981). 

In these flames, nitric oxide originates from the molecular 
nitrogen in air accord,-' ng to the Zeldovich mechanism (described 
elsewhere; e.g., Bowman, 1972). Some "prompt” NO formation is 
observed near the flame front at a few millimeters above the flame 
holder (considered as the main reaction zone). This is 
particularly noticed in fuel-rich cases where the Zeldovich 
mechanism cannot account for the measured NO^ levels. Several 
plausible reaction schemes for the prompt NO formation in rich 
flames involving the attack of hydrocarbon fragments on N 2 and the 
conversion of the products to NO have been suggested (Fenimore, 
1971; Bachmaier et al., 1973; Haynes; 1977b; Hayhurst and Vince, 
1983). In spray combustion, the presence of liquid droplets 
increases hydrocarbon cracking and thus intensifies the prompt NO 
production. 

The consumption of NO^ in the post flame regions appears to be 
independent of the fuel type or stoichiometry* Equivalence ratio 
seemed to only affect the NO* reduction rates. For 37 pm n-octam 
droplets, this accounted for a drop of 22, 27, and 31% at 0=0.85, 
1.0, and 1.2, respectively. Duterque et al . (1981) reported 
similar behavior in a stirred reactor. 

Changes in flame temperature were also found to affect the 
removal rate of N0^ in the post flame regions. Experiments with 
synthetic oxidizers were carried out at an equivalence ratio of 


1.0 (see section 4.5 for the experimental details). Figure 4.23 
compares the axial NO* and temperature profiles obtained for 37 pm 
n-heptane droplet combustion with various oxidizers. The data 
show that raising the flame temperature increases the NO* 
formation markedly , but also augments its downstream decay to 
constant levels. 

A possible explanation for the axial NO* reduction is the 
removal of NO by nitrogen atoms which are generated from the 
molecular nitrogen at high temperatures. Other reactions that can 
also contribute to the consumption of NO* are discussed in the 
following section. 

4.7.2 Nitrogen-Containing Fuels 

Several operating conditions were selected to study the NO* 
formation in spray combustion from fuel-N. Centerline NO, NO* and 
temperature measurements for stoichiometric combustion of 
n-heptane doped with 1% nitrogen by weight as pyridine, are shown 
in Figure 4.24. Effect of equivalence ratio for a fixed droplet 
diameter of 37 pm is demonstrated in Figure 4.25. Again, in the 
flame zone, NO and NO* rise quickly with temperature to their 
maximum values, and then fall slowly to a plateau in the post 
flame regions. Both NO and NO* peak faster at higher fuel-N 
contents and stoichiometries, indicating a strong dependence on 0 % 
concentration but a weak dependence on flame temperature for 
fuel-N transformation. Similar results for* different fuel-N 
concentrations are shown in Appendix B. 

Surprisingly, in each case, a lower percentage of the measured 



























peak NO was consumed far from the burner, as compared to the 

x 

results with nitrogen-free fuels. Table 4.2 summarizes the 

reduction percentages for the conditions studied. One may reason 

that the NO generated from the fuel-N is far greater than the 
x 

thermal NO and therefore its removal by the flame radicals is 
x 

less efficient. 

At this point a closer look at some of the processes that 

control the NO formation from fuel-nitrogen may help the 
x 

interpretation of the results. It is generally believed that, in 
fuel rich flames, the overall reaction mechanism for fuel-N 

conversion to NO occurs according to the following path 

x 

(Fenimore, 1976b; Haynes 1977a; Chen and Malte, 1984): 


Fuel-N -> HCN NH 



Among different reactions that can lead to the oxidation of the 

NH: species (N,NH,NH ,NH ) and NO formation, the reaction: 

1 2 3 

N+OH -» NO+H 


is considered to be the most important (Haynes, 1977a). The 
unoxidized NH^ species can react with NO as it forms to generate N 
by any of the following reactions: 


N+NO -» N +0 

2 

NH+NO -> N +0H 
2 

NH +N0 -* N +H 0 
2 2 2 . 

For lean mixtures, due to the abundance of 0 and the presence 

2 

of 0 atoms, the fuel-N transformation to NO is much easier. The 
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Table 4.2. Percent NO* consumed In the post flame gases of 
37 pm n-heptane droplets doped with pyridine. 


% N By Weight of Fuel 

(j)*«0.85 

4>-i .o 

<j>«1.2 

0.0 

19 

28 

20 

0.1 

9 

18 

20 

0.5 

5 

15 

20 

1.0 

4 

15 

20 


likely reactions are perhaps those suggested by Appleton and 
Heywood (1973), by which radicals such as NH and CN oxidise to NO 
according to: 

NH+O -4 NO+H 
NH+O OH+N 

N+O -4 NO-fO 

2 

and, 

CN+O -4 NO+C 
CN+O -4 CO+N 
N+0 2 -> NO+O, 

In ammonia oxidative systems, NO can also form when nitroxyl 

hydride (HNO) which is a product of amine species oxidation » is 

stripped from its hydrogen (Dasch and Bllnt,1984). This can 

happen as a consequence of radical termination reactions: 

HNO+NH^ -* NO+NHg 

HNO+O -4* NO+HO_ 

2 2 

or, the radical branching reaction; 

HNO+M -> NO+H+M. 

Other routes for NO formation and destruction in these flames 

involve OCN radicals (Crowhurst and Simmons, 1983, 1985a, 1985b) 

OCN+O CO+NO 

OCN+NO ^ N +C0 , 

2 2 


as well as 


GN+NO N +C0 
2 


CHAPTER 5 


Analytical Studies 


5.1 Introduction 

So far, It has been shown that, In spray combustion, 
variations in fuel properties, air preheating temperature and 
oxidizer constituents affect the NO^ production. These parameters 
are examined by an analytical model which is described in this 
chapter. 

The model is in two parts and consists of spray vaporization 
and spray combustion. First, droplet diameter and extent of 
vaporization before the flame are evaluated by the fuel spray 
vaporization model. Then, the calculated results are used as 
initial conditions in the spray combustion model that also 
considers droplet interactions and the Zeldovich mechanism of NO 
formation. Temperature and fuel, oxidizer and NO mass fractions 
around the burning droplets are determined from this analysis. 
Single component nitrogen free fuels are studied; therefore, 
atmospheric nitrogen is the only source of NO^ formation. 

Details of the assumptions and the governing equations as well 
as the comparison between the experimental and computed results 
are discussed in the following sections. Calculations of the 
tbermophysical properties are given in Appendix C. 


5.2 Spray Vaporization 

Given the initial state of a monodisperse fuel spray at the 
aerosol generator, this model predicts the final* droplet diameter 
and the extent of vaporization of the spray prior to the flame. 
Various equations describing the motion of droplets in the gas 
flow field, as well as energy and mass conservation equations, 
were considered. A stepwise droplet tracking technique was used 
to estimate the local vaporization parameters within the spray and 
to solve the coupled equations. 

S 

Starting with the equation of motion (Beer and Ghigler, 1983), 
the forces acting on a droplet travelling in a gas flow field are 
balanced according to 


<t *"Ss 3 jr =z ' /3 ”$ 's3 +c o 


*■ 0 

nr s p 

s 3 2. 


>2. 


(5-1) 


The droplet acceleration term on the left side of equation (5-1), 
which is determined from the rate of change of spray velocity due 
to divergence, is balanced by the gravity, buoyancy, and drag 
force terms on the right side. Considering a drag coefficient of 
C Q - 24/Re 0 (Re 0 <l), U', the relative velocity between the droplet 
and the surrounding gas can be determined from equation (5-1) 

u , 3<g-g>»jrH 

18 /x 


' (5-2) 
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Since equation (5-1) is one dimensional, its use in the model 
requires some justification. From visual observations of the 
spray pattern in the experimental setup (Figure 3.1), it was found 
that the droplets are contained within an inverted cone created by 
the dispersion jet. The spray cone is diluted by rapid dilution 
air entrainment near the baffle plate and diverges gradually to a 
maximum spread of about 1.5 cm along the 17.8 cm path to the flame 
holder. Because of the small cone angle, the radial variation of 
droplet size is expected to be negligible. Therefore, a one 
dimensional model should adequately describe the droplet motion in 
the preflame regions. Furthermore, the entrainment effect was 
determined experimentally (Cernansky and Sarv, 1983) and used to 
calculate the droplet travel time. This was done by integrating 
the ratio of the unit volume of the spray cone to the volumetric 

i 

flow rate of the combined dispersion and dilution air flows 
(described empirically as a function of axial distance). Also, 

the gas phase temperature in the model was adjusted based on the *■ , 

h 

dilution air entrainment results. 

At any given time, the total mass of fuel is constant. This 
can be expressed as 

m. + mp =* m + m, • (5-4) 

Assuming that the heat and mass transfer rates around the 
droplets are sufficiently fast to rapidly modify the local 
properties of the gas phase, the mass fraction of fuel vapor, Yp , 
is given by (Law, 1977) 

-(# 


where, Y c and Ya are the initial mass fractions of vapor and 

'o a 

liquid fuel in the combined vapor/gas system,, and 6 is the 
fraction of liquid vaporized defined as 



(5-6) 

(5-7) 


For a control volume of aerosol containing both liquid and 
vapor, the energy balance is given by 

5V < m 4-'5) L+( 'V”P c £5 v (n k-V c z <5-£> 

= "V„ ^ + t '”9 0 + ■ (5-8) 

Solving for from the above equation, the gas temperature can be 
expressed as 


"5 


= Vs ~k±k [hi 'll •‘Vs l . (5 _, 


c %( i+e \) 


(5-9) 


The evaporation rate for a droplet influenced by the relative 
gas phase motion is given by (Rao and Lefebvre, 1981) 


k, 


<7. sr 


m “ fr- & ('+6) {l + o.i-S'Re ) 


'5 


(5-10) 


where, and Cp are the conductivity and specific heat of the 

3 

gas surrounding a droplet of radius r^ , and B is the vaporization 
mass transfer number given by the following expressions 


6 = 


V 


(5-11) 


Assuming that the fuel mass fraction at the droplet surface is 
that corresponding to the saturated vapor pressure, one can 
integrate the Clausius-Clapeyron equation and use Dalton’s law to 
obtain 


‘ l 1 '^ a ' expt 7 <1/Ts ' 1/Tt)I ] 


m can also be expressed as 


(5-12) 


d t. 


(5-13) 


Substituting equation (5-3) into equation (5-10), and using the 
equality between equations (5-10) and (5-13) , we get 


3 /z . 


r s da_ Ci( „ c/s -) 


where 


and 


C. s O.l! 


18 £ J 


o. S" 


C -21 bn (1+8)' 

P a C P 

3 P 3 

The parameters in the equations given for Cj and are 


(5-14) 

(5-15) 

(5-16) 


temperature dependent. However, the gas phase and the droplet 
surface temperature are also related by 


W T i) 


(5-17) 


and T c is solved using the equality between equations (5-11) and 


(5-17) (Lewis Number “1). 

The droplet temperature must 
phase temperature calculations, 
droplet varies both radially and 
equation takes the form 

- r 1 ^ ) 


be determined prior to the gas 
Since the temperature inside a 
with time, the heat conduction 

(5-18) 


with initial conditions 



and boundary conditions 



(5-19) 


(5-20) 


t)=r (5-21) 

which has to be solved to get the temperature distribution within 
the droplet. A dimensionless radial coordinate, r=r/r s is used to 
include changes of droplet size during heating or cooling. 
Substitution of r into equations (5-18) to (5-21) yields 


yr 

Vmw< 
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with T\.(r,o)« T /(r) 

T t Q 

(5-23) 

(®I) 

(5-24) 

d r r=o 


and ■ T. . 

r* . S 

(5-25) 


The above differential equation can be easily replaced by a 
finite difference equation. In the explicit form as described by 
Law and Sirignano (1977), droplet temperature can be determined by 
simple marching. This technique requires small time steps to 
avoid instability. On the other hand, the implicit procedure is 
unconditionally stable and, therefore, is a more appropriate 
choice Using Ar and At to denote the radial and time intervals, 
the implicit finite difference approximation of equation (5-22) 
becomes 


T f u , 

(5-26) 

with T r , » Tr . , 

V V 

(5-27) 

T / * T /' » 
■*1 r z 

(5-28) 

where , 


A = 

(5-29) 

r * 

and 


. o( A fc 

V - — 2_* 

r/(Arf 

(5-30) 

The superscript prime specifies the new value of 

temperature 


determined at the end of a time step and the subscript i 


* » 


i: 


designates the radial index for N r meshes. Equation (5-26) 
represents a system of linear equations wich can be written in 
matrix form with ( £-># ), (1+2#) and (-^-#) becoming the 
coefficients of a tridiagonal matrix. Knowing the initial and 
boundary conditions as given by equations (5-27) and (5-28), the 
simultaneous system of equations is solved for N r -1 unknowns, TJ . 
Once the radial temperature distribution has been determined, the 
average droplet temperature can be evaluated by spatial 
integration. 

J 7> (in?)# 


T f « 


y.rui T 


(5-31) 


Substituting for from equation (3-31) into equation (5-9) gives 
the gas phase temperature, Tj . 

Equation (5-14) can now be integrated as a function of time, 
subject t:-' r„ (t»0) « r s . Dividing the total droplet travel time 
into time steps, the local temperatures, mass fractions, gas 
properties and droplet diameter can be evaluated at the end of 
each time step. The calculated values are then used as the 
initial conditions for the next step. This procedure is continued 
until the total travel time is reached. 


5.3 Spray Combustion 

It was noted that, in the present experimental setup, the fuel 
spray enters a self-supporting, one dimensional flame, stabilized 
on a screen flame holder. Due to vaporization along the spray 
path and droplet impactions on the screen wires, part of the fuel 


is prevaporized and premixed with the oxidizer. The 
prevaporization due to simple evaporation is calculated by the 
vaporization model just described, while the screen impaction 
contribution to prevaporization must be determined experimentally 
and entered manually. 

After the fuel spray enters the combustion zone, solution 
involves a transient-diffusive-convective problem. The basic 
formulation is similar to that of Botros et al. (1980) in their 
numerical treatment of a single fuel droplet burning in a reactive 
gaseous system. 

With assumption of uniform pressure, zero relative velocity, 
and one-step reaction kinetics for fuel oxidation, the governing 
equations in a spherically symmetric system reduce to 
continuity : 

( 5’ rrZ) ’° < 5 " 32 > 


fuel : 




oxidizer : 


X (P0r'^') ! .« u (5 . 34) 

9 St Kr' r 1 ~>r v 6 ' ° 0 ^ w> 


energy : 


pc p + (k r 1 w ah , 5 _ 351 

9 a r ' r l )r 3 l>r' F F ^ 3t>) 


state : 


P = 

where 

60 - 


/ 



(5-36) 

(5-37) 


is the mass diffusivity, u is the convective velocity due to 
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droplet vaporization and temperature variations in the gas phase, 

V is the stoichiometric coefficient, A is the pre-exponential 

It 

terra and E is the activation energy (for n-octane, A-4.6xl0 , 

12 

E»3C)C)00 and for methanol, A»3.2xl0 , E-30000 cal/mole; Westbrook 

and Dryer, 1984). 

Boundary conditions at the droplet surface (r»r ) are 

O 



(5-38) 


(5-39) 


(5-40) 


(5-41) 

r(.r,r t ),T s . 

(5-42) 

At the midpoint between the droplet centers 

t the ambience boundary 

conditions are 


(23) = (23) = (23) . 

N Dr '<» V Dr N Dr co 

(5-43) 

The initial conditions at t=0 are given as 

radial profiles 

Tg(r.O) - \(r) 

(5-44) 

Y F (r,0) = Y % (r) 

(5-45) 

Y o (r,0) - Y 0 (r) 

o 

(5-46) 

v(r,0) =* v Q (r) 

(5-47) 

with r s (0) ■ r s . 

(5-48) 



Applying the dimensionless coordinate, r, the governing 
equations ^5-32) to (5-35) as well as the boundary and initial 
conditions (Eqs. 5-38 to 5-48) can be represented as 


+ 2- 

^ r ' 

P fLJ i (p 2>r ^ ) = --VW, 

3^ r& 2 >r ' r r/ ~df 3 ' F 

p / Vfo , ¥_~ d ' 1 ' 0 )- _L_ (p Dr x & 

^ r 5 f r/ ' 00 


(5-32)' 


« (5-33) ’ 


(5-34) 


pc. (13+ ¥:Vj 

3 r 9 2t r s 




5-35)’ 




°aP_ ( 

r s 


(Ht) + P q (ir-'iH) ( /- Y ) = o 

-)r 'f 9 s it 1 F < 

- P (ir t£)Y 0 -o 

^ t ?a ) -p (v- h)L = o 
If 3. s dt J 


c> r r r I 


n. = t 

r= I 5 




~aV c 


( a ’p \ _ / d t o n / yjs \ 

tr V s <p " ^=oo" k ^ ; r= 

T 3 (.r>°) = T 3o (t) 

f (r,o) = Yp tr) 


00 


(5-38)’ 

(5-39)’ 

(5-40)’ 

(5-41)’ 

(5-42)’ 

(5-43)’ 

(5-44)’ 

(5-45)’ 
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V 0 (r, 0 ) = Y 0 o (^o) 
v(.r,o) = tr o (.r,o) 

r s ( o) --■ % . 


(5-46)' 

(5-47)' 

(5-48)' 


A finite difference approximation scheme was used to solve 

these governing partial differential equations. For better 

resolution in the high gradient regions around the droplet 

surface, the radial spacing was ramped according to 

*i * *i-i +0 '2+0.06(1-2), i>2. (5-49) 

Also, in order to apply this model to sprays, the burning rate was 

corrected to account for droplet interactions. Since the droplet 

d t 

surface regression rate, XT, - •— is directly proportional to the 

6 dis 

mass burning rate, such correction was accommodated by adjusting 
the velocity term in the model. Calculations of the correction 
factor are discussed in the next section. 


5.4 Droplet Interactions 

A mathematical procedure haa been developed to transform the 
quasi-steady transport equations describing the Stephan flow field 
around the interacting droplets to the Laplace equation. Details 
of the analysis as well as the underlying assumptions have been 
given by Labowsky (1976-1980) and therefore, will not be repeated 
here. However, the following discussion should provide some 
insight to the development of the model. 

The combined fuel/oxidizer conservation equation can be 
written as 



(5-50) 


Defining the velocity potential, ^,as 


v c w c 




- ftl L ( '"'’f + ( o’^W 0 ) / < l+Y o,® U3 


'O w O 


JZt, ( 1 + 8) 


(5-51) 


the mass averaged gas velocity, v, can be expressed according to 
V - JDtn (HB)V^ (5-52) 

where B is the combustion transfer number given by 


B 


y f +y 

' re n 


S O/P 


~%>^o 



(5-53) 


Equation (5-52) together with the mass conservation equation 
(V. Pt/'-O) yields the Laplace equation 

ty -0. (5-54) 

From the definition of (Eq. 5-51), the boundary conditions 
become 

r 1 on the surface of droplets 

v H 

10 far away from the droplets. 

The above equation is solved using a superposition technique 
called the method of images. Briefly, in order to solve for the 
fields, a series of point sources and sinks are placed inside 
every droplet in the array which satisfy the boundary conditions. 

The burning rate of the droplet k in an array is determined 


by integrating the mass flux over its surface area 


Substitution for v from equation (5-52), gives 


i -4ncp®in (1+B) — ! — IvU'.dS. 
k 5 3 Hnr, " 


(5-56) 


The right hand side of equation (5-56) can be rewritten as 

K ■ • (5 - 57) 


where is the burning rate of droplet k if it were burning as 
an Isolated droplet and ^ is the burning rate correction factor 
for droplet interactions defined as 


7. 


H nr* 


V W .dS, 


(5-58) 


For an array of N interacting droplets, 7 can also be shown 


as (Labowsky, 1980b) 


7-Z « kt Vi 

'k f.-i ' 


(5-59) 


Where o( designates elements of an NxN matrix representing the 
transfer rates of droplet k, when it is in the field of a unity 
potential droplet, £ , with the rest of the array being at zero 
potential. Using the method of images, these elements can be 
determined by adding the images source strengths inside the k-th 
particle that are associated with the $ field. 

As it appears, the calculations of the combustion rate ' 
correction factor are independent of the fuel or oxidizer 
properties and only depend on array configuration. 


5.5 Array Geometry 


In this work, the correction factor is computed for an evenly 
dispersed configuration Inside the flame holder. The droplets are 
divided into geometrically similar planar arrays* separated by the 
Interdroplet spacing, as shown in Figure 5.1. Except for the 
droplets near the burner wall, every droplet is encircled by six 
immediate neighbors placed at the corners of a hexagon. For a 
burner diameter, d, and center-to-center droplet spacing, JL , the 
number of droplets that can be placed on a layer is equal to 


number of droplets per layer 


71 d ? '/4 

C/z Si " 60 


(5-60) 


The total number of droplets per unit volume or number density 
is given by 


droplet number density< 


O^Zr) ( Sin 60 ) 


71 


(5-61) 


Since the burner dimensions are large compared to the droplet 
spacings (i.e. h/£»l), the above expression simplifies to 

i 

Z 

droplet number density « — - (5-62) 

<£ sin 60 

Usually, droplet number density is measured by optical 1 
counters or determined from the operating conditions (as in 
well-chare Ktarized sprays). Thus, droplet spacing in a uniformly 
distributed raonodisperse spray, can be calculated from equation 



77 


(5-61). 


5.6 NO Kinetics 

Nitric oxide formation reactions involving the oxidation of 
atmospheric nitrogen (Zeldovich mechanism) are considered. These 
reactions and their rate data are listed below (Bracco, 1973; 
Bowman, 1975; Peterson and Laurendeau, 1982; Dasch and Blint, 


1984). 


0 +N 

2 

N +0 


0 +N 
2 2 


Forward rate constant 

, 6 


Backward rate constant 

O 


NO+N 1.36xlo"Up(-75,400/RT) 3.1xl0 13 exp(-334/RT) 

9 o 9 o 

N0+0 6.4x10 T exp(-'6260/RT) 1.5x10 T exp(-38,640/RT) 


0+0+N 6.0xl0 3 (118,000/tfT) 
2 

x exp(-l 18 ,000/$T) 


14 


1.0x10 


The rate constants for the bimolecular and termolecular reactions 

3 6 2 

are in units of cm /mole sec and cm /mole sec, respectively. In 

view of the fact that NO transformation to NO is not considered, 

2 

these reactions account for the total NO production. 

x 

Assuming steady-state N atom concentration and equilibrium 

concentration for 0 atom with respect to 0 , concentrations of N 

2 2 

and NO are determined as their species conservation equations are 
solved simultaneously with the energy, fuel and oxidizer 
equations. The initial and boundary conditions are similar to 
those for oxygen as discussed in section 5.3. 


5,7 Representative Results 

5.7.1 Preflame Spray Vaporization 

Calculated results for stoichiometric sprays of methanol and 

C - 3 - 


n-octane are presented in Figures 5.2 and 5.3 for an initial 
droplet diameter of 50 pm and droplet and gas temperatures of 24 
*C. These plots show changes in droplet diameter, gas phase and 
droplet, temperatures, extent of prevaporization and the mass 
transfer number during vaporization. The last point on each curve 
corresponds to conditions at the end of the droplet travel time in 
the preflarae regions of the combustor. 

After vaporization is initiated, the droplet surface 
temperature drops to its lowest value. Then, as illustrated in 
Figure 5.4, it increases gradually and approaches the falling gas 
and droplet temperatures. The reduced gas temperature and 
increased fuel vapor concentration decrease the rate of 
vaporization substantially, Inhibiting further droplet size 
reduction. This is also demonstrated in the plots of the mass 
transfer number versus droplet travel time, where B number 
decreases continuously and nears zero as the fuel vapor 
concentration reaches its saturated value. 

When vaporization occurs adiabatically, heat is extracted from 
the liquid fuel as well as the gas phase to vaporize the droplets. 
Among the fuels used in this study, methanol and isopropanol 
require higher heat of vaporization than n-heptane and n-octane. 

At the same time, the air to fuel mass ratio for a fixed 
stoichiometry is less for the alcohols as compared to the alkanes, 
leading to a more effective cooling and enrichment of the gas 
phase by droplet vaporization. Under such conditions, the vapor 
pressure at the surface of the fuel droplets decreases and thus 
hinders further changes in the evaporation characteristics of the 
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fuel spray, even for long droplet residence times. 

In these calculations, the spray is considered to vaporize 
adiabatically while travelling inside the upper and lower housing 
assemblies. But, as the fuel aerosol enters the flow reducing 
sections, additional heat is transferred to the cold spray from 
the combustor walls which are kept at the air preheat temperature. 
Heating of the aerosol by constant temperature walls is treated as 
a thermal-entry length problem. This results in an increase in 
the gas and droplet temperatures and a decrease in droplet size as 
shown in Figui-es 5.2 and 5.3. In the absence of the wall heating, 
adiabatic spray vaporization continues for the entire droplet 
travel time ms shown by the dotted lines. 

Table 5.1 summarizes the calculated vaporization parameters 
for the fuel sprays tested in this study, based on an initial 
droplet diameter of 50 pm and different air preheating 
temperatures. As expected, the final droplet size and the extent 
of prevaporization are strongly influenced by fuel volatility and 
ambient temperature. These predictions are used later as initial 
conditions in the spray combustion model. 

Figure 5.5 compares the predicted results of the aerosol 
temperature with the actual measurements made along the spray 
centerline. Since droplet impaction and agglomeration on the 
thermocouple tip was unavoidable, it is assumed that the 
measurements represent the average of the droplet, gas phase and 
droplet surface (wet bulb) temperatures. Both measured and 
calculated values correspond to initial droplet diameters equal to 
the N0 X minima, air temperature at 22 *C and overall stoichiometry 
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of 1.0. In spite of the uncertainties in the measurements and the 
model, predictions, the trends are similar. 

5.7.2 Burning Rate Correction Factor 

Effects of droplet number density and interdroplet spacings on 
spray burning rate are examined for evenly dispersed arrays of 
raonosized droplets (similar to the one shown in Figure 5.1). In 
such configuration, droplets that are located around the array 
center experience a stronger interaction than those located close 
to the outer edge. 

Computed burning rate correction factors for three arrays of 
interacting droplets are listed in Table 5.2. The calculations 
were done for a fixed interdroplet distance of 20 droplet 
diameters (<£**40). In all cases, the burning rate correction 
factor, y , for any given droplet is less than unity, indicating 
slower burning compared to that of an isolated droplet of the same 
size. 7 increases as one moves away from the center of the array 
and decreases with increasing number of droplets,, 

For closely packed arrays or those that contain large number 
of droplets (>37), the accurate determination of 7 b 7 Labowsky’s 
droplet superposition procedure becomes exceedingly tedious. 

Thus, A second superposition method was used to extend his 
technique to cases that involve a great number of droplets. This 
development is based on the fact that any monodisperse array can 
be constructed by superimposing arrays of six droplets around a 
primary array of seven. For instance, a 55 droplet array can be 
composed of eight consecutive arrays spaced about the center of 
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TABLE 5.2 Computed burning rate correction factors for three 
monodisperse arrays of droplets with £ -40. 


Number of Droplets 
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Array Average 
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the primary array* aa demonstrated in Figure 5.6. With every 
addition, the burning rate of the entire configuration reduces by 
a factor which is equivalent to the burning rate correction factor 
of the superimposed array. Figure 5.7 shows the calculated 
average burning rate correction factor as a function of droplet 
spacing for various raonodisperse arrays. The data points and the 
solid lines represent the results obtained from the droplet and 
array superposition techniques, respectively. The results are in 
good agreement with a maximum difference of 1.0% f >r the studied 
cases. 

When droplet interactions are considered in a spray combustion 
model, the burning rate correction fa '.tor has to be reevaluated as 
the interdroplet spacing changes. However, successive application 
of the droplet interaction model could add up co a considerable 
increase in the computation time. To overcome this problem, the 
dependence of the burning rate correction factor on £> is expressed 
as a functional relationship. Computed values of 7 are used to 
determine the ’best fit’ in the form of four degree polynomials. 
These expressions are utilized in the spray combustion model to 
facilitate the evaluation of 7* 

5.7.3 Spray combustion and NO formation 

Computations were performed for stoichiometric sprays of 
n-octane and methanol. In all cases, the initial conditions at 
t-0 were determined by the spray vaporization model. The time 
domain was set equal to the duration, after which, the calculated 
results were at their steady state values. 







To initiate combustion, the gas phase temperature was raised 
to 1100 K by reacting a part of the fuel vapor and the oxidizer. 
Further fuel and oxidizer consumption as well as the temperature 
rise due to the heat release are determined by the spray 
combustion model. The continous increase in the gas phase 
temperature eventually leads to droplet ignition which is seen as 
a bulge in the temperature profiles. Figure 5,8 to 5.10 show the 
development of temperature and species profiles for 27.0, 33.2 and 
47.7 pm n-*octane droplets, corresponding to initial droplet 
diameters of 37.0, 43.3 and 57.5 pm, respectively. The time 
dependence of the profiles are noted on the plots. Followed by a 
rapid reduction in the fuel vapor concentration, the flame 
temperature reaches its maximum value close to the droplet 
surface. The location of the maximum temperature coincides with 
the crossing of the fuel and oxidizer profiles and moves with 
changes in droplet vaporization rate and oxygen concentration 
during combustion. For small droplets, gas phase combustion, and 
droplet ignition are intensified by the high preflame extent of 
vaporization. Much of the NO is formed in the hot gases 
surrounding the burning fuel droplets, where preraixed type of 
burning appears to dominate. In the case of large droplets, 
temperature rises slowly as the reactions occur in a relatively 
narrow zone. This is attributed to a lower droplet vaporization 
before the flame front and the slow propagation of the fuel vapor 
emerging from the surface of the burning droplets. After the 
initial temperature rise, diffusional burning is also observed 
near the droplet surface, as indicated by the peak in temperature 
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Figure 5.10. Evolution of temperature and mass fraction profiles 
for a 47.7 pm droplet in a raonodisperse spray of 
n-octane at (J)= 1.0 and £ = 0.43. 




profiles. However, NO production for large droplets increases as 
they achieve elevated flame temperatures due to their higher 
Damkohler number. At some intermediate size, where minimum NO* is 
predicted, the flame temperature and the fuel vapor concentration 
in the gas phase fall in between the values observed for the small 
and large droplets. Under this situation, droplet Interactions 
reduce the burning rate and lower NO formation by confining its 
production to a thin reaction zone. 

In these plots, there is little change in the nitrogen 
profiles as a very small amount of it reacts to form nitric oxide. 
A slight but gradual decline in the N 2 , levels is noticed as the 
product gases move out radially and change the nitrogen mass 
fractions. Since, nitrogen is always in abundance, nitric oxide 
formation is primarily influenced by the gas phase temperature and 
oxygen concentration. The location of maximum NO formation is in 
the vicinity of the peak temperature, but it moves outward to 
regions of higher oxygen concentration as the temperature peak 
begins to diminish. 

Similar profiles were also obtained for methanol sprays. The 
main difference between the methanol and n-octane results is in 
the relative location and the magnitude of the flame temperature. 
For methanol, the flame front stays very close to the droplet 
surface due to its high heat of vaporization. This phenomenon has 
also been recognized by Kesten (1972). Relative to n-octane, 
methanol has a lower heat of combustion which gives a lower flame 
temperature. These two events result in a lower NO formation for 
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methanol. In addition, inclusion of the prompt NO formation 
mechanism can te expected to increase the difference, by 
predicting a higher contribution for n-octane. 

Predicted and experimental results for n-octdne, showing the 
effects of droplet size and air preheating temperature on NOy 
formation, are compared in Figures 5.11. Although the calculated 
results are higher than the measured values, the general trends 
are In reasonable agreement. As the air preheat temperature 
increases, more fuel is vaporized before the flame front, leading 
to a rapid increase in the flame temperature and the weakening of 
the droplet interactions. The model also predicts the shifts of 
the NO^ minimum point towards larger droplet diameters with 
increased air preheating temperatures. In these analyses, the 
contribution to the total extent of preflarae vaporization due to 
screen impaction was not considered. Such consideration requires 
experimental measurements of fuel droplet losses. Inclusion of 
this effect is expected to reduce the skewness of the predicted NOy 
trends on both sides of the minimum NOy point, thereby improving 
the agreement between the model and experiment. 

In summary, the computed results indicate the importance of 
the extent of prevaporization and droplet interactions on NOy 
production. Gas phase stoichiometry has a significant effect on 
the volume of the hot gases surrounding a fuel droplet, where NO 
is formed. On the other hand, the release of fuel vapor from the 
droplet to the reaction zone is controlled by droplet interactions 
which reduce the volume of the hot gases and NO production by 
bringing the flame closer to the droplet surface. 
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While a qualitative agreement exists between the measurements 
and the predictions, the model overestimates the NO enijsion 
levels. Higher NO predictions are directly related to the greater 
flame temperatures computed in the model, using a single step fuel 
oxidation rate. The discrepancies can be remedied by the 
following considerations. Use of detailed kinetic mechanisms for 
the fuel oxidation and NO^ formation can Improve the predictions 
significantly. An accurate determination of the actual droplet 
travel time in the preflame regions of the combustor is also 
desirable due to its importance in the evaluation of the initial 
conditions for this model. Also, the applicability of the 
spherico-symraetric governing equations if, restricted to zero 
gravity or convection free situations. Convective effects distort 
the spherical flame shape which can bring about different droplet 
interactions and NO^ formation. A detailed analysis of the flow 
field around the burning droplets would be necessary in order to 
relax the restrictions imposed upon spherical flames. 

These refinements should alleviate the differences and make 
the model a useful tool for defining the operating conditions at 
which N0 X formation is minimum. 


CHAPTER 6 


Conclusions and Recommendations 


6.1 Conclusions 

A one dimensional monodisperse fuel spray combustion facility 
was developed and used to study the formation of oxides of 
nitrogen under various operating conditions. Inflame and 
postflame measurements of temperature, NO and NO^ , and other 
emissions were made over a 26-74 pm droplet size range and 
different stoichiometries. Single component fuels including 
methanol, isopropanol, n-heptane and n-octane, as well as 
nitrogen-containing fuels were tested. 

Postflame measurements of the exhaust gases from the 
combustion of nitrogen-free fuels indicated the existance of an 
optimum droplet diameter for lower N0 X production in the 43-58 pm 
droplet size range. Up to 30% reduction in the N0 X emissions, 
relative to the small droplet limit, was detected. The observed 
behavior was attributed to the effects of droplet size, gas phase 
stoichiometry and droplet interactions in the flame. 

Extent of prevaporization determines the gas phase 
stoichiometry and affects the temperature and the volume of the 
hot gases surrounding the fuel droplets. The thickness of the 
reaction zone, however, is influenced further by droplet 
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interactions which control the s^ray burning rate. All of these 
factors depend on the initial droplet diameter and the fuel 
evaporation characteristics. For small droplets, rapid 
vaporization before the flame front intensifies gas phase 
reactions, leading to a fast Increase of temperature and NO^ 
production > Large droplets evaporate less, but form more NO as 

X 

they attain higher flame temperatures during their lifetime. 

Droplet interactions are weak for both small and large 

droplets. This is caused by rapid preflame vaporization and 

increased droplet spacings in the former case, and fewer 

interacting droplets in the case of large droplets * The optimum 

droplet diameter for minimum Nf^ formation occurs at some 

intermediate size with relatively stronger Interactions. 

A mathematical model for the analysis of preflame droplet 

vaporization and NO formation in the combustion of nitrogen-free 

x 

fuel sprays, was developed. The predicted results indicated a 
qualitative agreement with the measurements and demonstrated the 
importance of the evaporation and combustion characteristics on NO 

X 

emissions. 

With nitrogen-containing additives, the nitrogen 
transformation process to NO^ was less dependent on flame 
temperature and was influenced primarily by oxygen concentration. 
At high fuel-N contents, NO decreased continuously with 

X 

increasing droplet size. For small droplets, oxidation of the 
pyrolysis fragments of the nitrogen additive improved as a result 
of higher prevaporization and better oxygen penetration through 
the droplet flame zone. 
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6.2 Rec ommendations for F uture Work 

■W - ■' — — ■ ■. .i , .— , l—i ■■■■■■■.,.! i. ,i i ■ i - —I— — 

While this study points out the potential for minimizing NO 

x 

emissions from fuel spray combustors, further experimental and 
theoretical advancements are still needed before 'the acquired 
knowledge can be applied to ouch systems. 

Suggested experimental studies and their objectives are listed 
below. 

(a) In situ droplet sizing and number density measurements in 
the flame and prior to the flame is required to obtain 
valuable information on vaporization and burning rates of 
interacting fuel droplets. The number of interacting 
droplets that survive the screen impaction and the total 
extent of preflame vaporization can be determined from 
such measurements. 

(b) Detailed analysis of combustion products and probe 

reactions are necessary to explain the observed NO 

x 

consumption in the post flame regions. 

(c) Polydlspersed sprays and the effect of Sauter Mean 

Diameter on NO formation should be considered for their 
x 

applicability to practical spray combustors. 

(d) Inflame and postflame analyses of the nitrogenous species 
must be carried out to get a better understanding of 

nitrogen transformation and NO yield in fuel sprays that 

x 

contain chemically bound nitrogen. 

As far as the modeling aspect of spray combustion and NO 

x 

formation is concerned, improvements can be achieved in several 
ways. 


(e) A more detailed chemical model for fuel oxidation i$ 
needed to improve prediction of NO production. 

(f) Flame shape distortion due to convective affects and its 

impact on droplet vaporization, droplet interactions and 

NO formation requires some consideration. An in-depth 
x 

analysis of the flow field around the burning fuel 
droplets is needed to account for such phenomena. 

(g) Fuel-N conversion to NO should be studied by implementing 

x 

a comprehensive species transformation mechanism into the 
spray combustion model. 
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APPENDIX A 


Supporting Experimental Facility 
A.l Introduction 

The major components of the experimental setup, i.e., the 
monodisperse aerosol generator and the temperature and NO^ 
measurement systems, were described in Chapter 3„ This appendix 
provides a description of additional facilities including a fuel 
prevaporizer and other exhaust gas analyzers. 

A. 2 Fuel Prevaporizing System 

In order to determine the small droplet size limit, a fuel/air 
heating arrangement was added to the experimental facility. This 
set up provided prevapor ized/preraixed fuel/air mixtures by 
combining the liquid fuel and air in a tee and heating the mixture 
in a stainless steel tube. The tube was formed into a vertical 
three-loop helix and wrapped with heating tape. The fuel vapor 
was heated to abc it 65 °C before entering the burner through the 
normal dispersion air path. This heating was sufficient to keep 
the mixture above the dew point of the fuel and to avoid 
condensation of the fuel vapor on the combustor walls. Normal 
cooling of the fuel vapor in the combustor resulted in a mixture 
temperature of about 26 °C at the flame holder exit. 


A. 3 Sampling and Measurements of Exhaust Gases (UHC, CO, CO 2 and 0 2 ) 
Inflame and postflame combustion gases were sampled and 
analyzed for hydrocarbons, carbon monoxide, carbon dioxide and 
oxygen concentrations. Hydrocarbons were measured as propane 
using a heated flame ionization detector (Beckman Model 402). CO 
and C0 2 were measured using nondispersive infrared analyzers 
(Horiba Models AIA-21-AS and AIA-21, respectively). 0 2 
concentration was measured using a polarographic analyzer 
(Teledyne Model 326AL). Due to the high flow requirements by the 
analyzers, the combustion gases were sampled under atmospheric 
pressure using a 3.17 mm (O.D.) stainless steel tube with a 2 mm 
opening at its tip. 


APPENDIX B 


Complementary Results 


B.I Introduction 

This appendix reports experimental results taken in addition 
to those presented in Chapter 4. In particular, the effects of 
fuel nitrogen content and droplet size on the NOy flame structure 
are demonstrated in section B.2, while the results of the exhaust 
emission measurements for a number of operating conditions are 
given in section B.3. 

B»2 Axial NO, NO^ and Temperature Profiles for Fuel-N Additives 

Representative results for the axial variations of 
temperature, NO and NO^ concentrations for n-heptane droplets 
doped with 1% nitrogen by weight as pyridine, were shown in 
Figures 4.24 and 4,25. Additional results at different fuel-N 
contents are shown in Figures B.1-B.5. After a rapid increase 
with temperature in the flame zone, NO and NOy reduce gradually in 
the post flame regions to their steady state levels. Results for 
different size droplets indicate that more fuel-N converts to NO* 
at smaller diameters. These general trends are similar to those 
discussed previously in sections 4.6 and 4.7.2. 
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Figure B.5. Centerline NO, NO* and temperature profiles for 37 pm 
n-heptane droplets doped with pyridine (0.5% N by 
weight). cj,b,x -<$)«« 0.85; A ,A,+ - 4) » 1.2. 







B.3 Exhaust Emission Measurements 


Exhaust gas analyses along the burner centerline were carried 
out In order to determine the local equivalence ratios and 
combustion efficiencies. Concentrations of major products and 
reactants Including unburned hydrocarbons, CO, C(^ and (J> were 
determined from the combustion products of n-octane aerosols for 
different stoichiometries and droplet sizes. The hydrocarbon 
sampling line was heated and therefore, the measurements are on a 

wet basis. The CO, CO and 0 measurements, however, are on dry 

2 2 

basis, since water was condensed and removed from the sampled 

gases before they were analyzed. 

Figure B.6 and B,7 show the effect of equivalence ratio on the 

axial variations of UHC, CO, CO^ and 0^ concentrations In a 37 j,un 

aerosol of n-octane. The results demonstrate a rapid increase in 

% 

CO and CO with a substantial reduction In 0 and UHC 
2 2 
concentrations in the flame zone. Further consumption of 0^, UHC 

and CO occurred in the post flame regions with CQ, values reaching 

steady levels at about 5 cm above the flame holder. For lean and 

stoichiometric mixtures, the exit plane CO measurements were very 

low (order of ppm’s), but increased for rich mixtures. On the 

other hand, the hydrocarbon emissions were consistently low for 

different equivalence ratios, indicating complete combustion for 

all cases. 


* The UHC emission measurements are shown only for the post flame 
regions. Large fluctuations in the UHC readings were observed 
near the flame zone where fuel droplets were drawn into the 
sampling probe, due to the high sampling flow requirements by the 
hydrocarbon analyzer. 
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. Effect of equivalence ratio on centerline CO, CO 2 , O 2 
and UHC profiles for 37 pm n-octane droplets. 
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The effect of droplet diameter on exhaust emissions was 
examined at an equivalence ratio of 1.0. Figure B.8 shows the 
UHC, CO, CO 2 and O 2 measurements at 10 cm above the flame holder. 
The highest CO and 0 2 , and lowest CO 2 concentrations were 
observed around 43 pm droplet diameter (i.e., the observed NO^ 
minimutn point for n-octane at (J)»1.0). This behavior is caused as 
a result of Increased droplet Interactions and suppressed flame 
temperature which lead to reduced 0 2 consumption and higher CO 
formation. Hydrocarbon concentrations appear to increase with 
increasing droplet diameter. 

Computed equivalence ratios, based on the exhaust gas 
analysis, are about 12% leaner than the measured input values. 
However, gas chromatographic analyses of the exhaust gases for CO 2 
and CO concentrations, revealed only a 3% difference in the 
overall stoichiometry relative to the measured input values, which 
is well within the accuracy of the measurments. The discrepancy 
in the measurements by the exhaust analyzer is due to 
instrumentation and calibration errors. However, this does not 
bias the above discussion, since, it is based solely on the 
relative changes of, concentrations rather than the absolute 


values. 


Figure B. 
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3. Droplet size effect on postflame CO, CO^, 0^ and UHC, 
for n-octane at (f> «■ 1.0. 
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APPENDIX C 


Thermophysical Property Calculations 


Certain properties of liquids, such as latent heat of 
vaporization and vapor pressure, are influenced primarily by 
temperature. For a mixture of ideal gases, transport properties 
vary not only with temperature, but also with concentration of the 
individual components of the mixture. These variations were 
included in the analytical computations. 

The temperature dependence of each property was established as 
a polynomial expression which fit the tabulated data with best 
accuracy. At any given temperature, for known species 
concentrations, the multicomponent transport properties are 
calculated based on the following formulae (Perry and Chilton, 
1973) 


density: 


P = 


specific heat: 



p 

R°T 2. Yj/vlj 


N 


Y: c f 


j=i J 


'/3 


thermal conductivity: 

w' 


absolute viscosity: 




- J* 


N , Vz 

2. x • M; W; 

J ' J J 


N' '/z 


These values were updated during the computations to account for 


changes in temperature and species concentrations. 


APPENDIX D 


Temperature Corrections 

As pointed out in Chapter 3, a bare Pt/Pt-13% Rh thermocouple 
was used for temperature measurements. The results presented in 
Chapter 4 indicate that all temperature readings and, in 
particular, the maximum flame temperatures are lower than 
expected. Several factors including convective heat losses to the 
gas sampling probe, conductive heat losses along the thermocouple 
wires, radiative heat losses from the thermocouple junction and 
catalytic reactions on the surface of the junction are considered 
to be responsible for the discrepancy. The latter case was found 
to be important only with new thermocouples. Catalytic activity 
on the thermocouple junction was tested in a manner similar to 
Kent (1970). After inserting a thermocouple into a flame, the 
flame was blown out, leaving the thermocouple exposed to the cold 
fuel/air mixture. It was noticed that while aged thermocouples 
cooled down quickly to the gas temperature, new ones continued to 
glow at a temperature of about 800 °C, indicating surface reactions 
on the junction. Reduced catalytic effects were observed after a 
few hours of operation at high temperatures. 

The actual temperature is often calculated from an energy 
balance on the thermocouple junction (Fristrom and Westenberg, 
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1965). This requires a knowledge of the local gas properties as 
well as the eraissivity of the junction, which are generally 
difficult to evaluate. As a result, any correction determined by 
this method is subject to a high level of uncertainty. 

Instead, a more direct approach was taken to estimate the true 
temperature. Four Pt/Pt-13% Rh thermocouple wires with different 
wire diameters (0.05 to 0.25 mm), but equal lengths (30 cm), were 
used to measure the flame temperature of a Bunsen burner at a 
fixed location. The bead diameter of each thermocouple was 
approximately 50% larger than the wire diameter. Figure D.l shows 
the variations of the flame temperature with the wire diameter. 

As expected, the losses are lower with smaller wires. In fact, 
the small heat dissipation resulted in the melting of the 0.05 mm 
thermocouple wires shortly after a maximum temperature of 1781 °C 
was measured. From these results, the actual flame temperature 
can be estimated by extrapolating the experimental trend to a 
thermocouple bead of zero diameter. With the 0.075 mm 
thermocouple, the combined conductive and radiative corrections 
are estimated to be about 72 °C at a measured flame temperature of 
1772 °C. 

Another source of error in the reported data is related to the 
interference of the stainless steel gas sampling probe. The 
probe, as shown previously in Figure 3.2, was adjacent to the 
thermocouple and is believed to have acted as a heat sink, . 
reducing the convective heat transfer to the thermocouple bead. 
Without the probe, the Bunsen burner flame temperature was about 
222 °C higher than that measured in its presence. 
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In summary, the total corrections for a measured temperature 
of 1550 °C could be as high as 294 °C. This amounts to an absolute 
temperature of about 2117*K which is close to the maximum 
adiabatic flame temperature. However, these corrections are not 
constant and require reevaluation at other temperatures. With 
respect to the volume of the temperature data collected in this 
study, determination of the actual temperature for all 
measurements was found to be impractical and, therefore, was not 
pursued . 


APPENDIX E 


Computer Programs 


E.l Droplet Interactions 


FORTRAN Variable Def Inltlon 


A (N,N) 

AL1 

ETA 

I 

IG 

I GEN 
J 

JJ 

K 

K1 

L 

N 

R (N) 

S ( IG,K) 

X < IG,K> 

XO (N) 

Y ( IG,K) 

YO (N) 

Z ( IG,K) 

ZO (N) 


Interaction matrix. 

Distance between source and particle I, 

Average array burning rate correction factor. 
Particle index. 

Index for number of source generation. 

Total number of source generation. 

Particle index. 

Particle index. 

Index for total number of sources per generation. 
Source index in previous generation. 

Particle index. 

Number of particles. 

Radius of particle I. 

Strength of source S. 
x-coordinate of source S. 
x-coordinate of particle I. 
y-coordir.ate of source S. 
y-coordinate of particle I. 
z-coordinate of source S. 
z-coordinate of particle I. 
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PROGRAM INTERACT 
C 

C AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

C A PROGRAM TO ACCOUNT FOR DROPLET INTERACTIONS IN ARRAYS 

C A OF MONODISPERSE OR POLYDISPERSE DROPLETS 

C A A A A AAA A A A A AAAA A A A AAA AAA A A A A A A A A A A AAAA A A A A A A A A A A A A AAA A A A A 

c 

DIMENSION X( 19 , 9000 ) ,Y(19,9000) ,Z(19,9000) ,X0<19> ,Y0(19) 
DIMENSION Z0( 19 ) ,A( 19,19) ,R(19) ,S( 19,9000) 

READ( 7 , A ) N 
READ <7,A) .IGEN 
DO 10 1=1, N 

10 READ ( 7 ,*) X0( I ) ,Y0( I ) ,Z0( I ) ,R( I ) 

DO 11 L=1 ,N 
A(L,L) =R(L) 

I=L+1 

DO 12 J *1 ,N-1 
X( 1 ,J) =X0(L) 

Y(1,J)=Y0(L) 

Z ( 1 , J) =Z0 ( L ) 

S( 1,J) =R( L) 

12 CONTINUE 

DO 13 IG = 2 , IGEN 
K1=(N-1)AA( IG-2) 

DO 14 K=1,(N-1)AA(ig-1) 

IF (I.LE.N) GOTO 50 
1 = 1 

50 AL1=SQRT( (X( IG-1 ,K1 ) -X0 ( I ) ) AA2+( Y< IG-1 ,K1 ) -Y0 ( I ) )AA2 
A +( Z( IG-1 ,K1 ) ~Z0( I ) > AA2 ) 

S(IG,K)=-R(I)AS( IG-1 ,K1 ) /AL1 
IF (I.NE.L) GOTO 60 
A ( I , I ) = A ( I , I ) +S ( IG , K ) 

GO TO 70 

60 A(L,I)=A(L,I)+S(IG,K) 

70 X( IG,K) =X0( I ) -R( I ) AR( I ) A ( (X0( I ) -X( IG-1 ,K1 ) )/ALl*A2) 

Y( IG,K) =Y0 ( I ) -R( I ) AR( I ) A ( ( Y0( I ) -Y( IG-1 ,K1 ) ) /ALIA A2 ), 
Z(IG,K)=ZO(I)-R(I)AR(I)A( (Z0(I)-Z(IG-1,K1) )/ALlAA2) 

IF (FL0AT(K)/(N-1) .NE.INT(K/(N-1) ) ) GOTO 80 
K1 =K1-1 
80 I=IA1 
14 CONTINUE 

13 CONTINUE 
11 CONTINUE 

WRITE (5,350) 

350 FORMAT ( ' THE INTERACTION MATRIX , ALPHA : ' ) 

ETA=0. 0 
DO 19 1=1, N 

WRITE (5,300) ( A( I , JJ) , JJ = 1 ,N) 

DO 19 J=1,N 
ETA=ETA+A( I,J) 

19 CONTINUE 


< < 


300 FORMAT (19F7.3) 

WRITE (5,400) ETA/N 

400 FORMAT </,' ARRAY AVERAGE CORRECTION FACTOR 
STOP 
END 
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E.2 Sprav Vaporization 


FORTRAN Variable Definition 


A 

(N) 

3 - Y. 

A1 


Locally defined variable. 

ADIAB 


Adiabatic /non adiabatic vaporization flag. 

AIRFLO 


Air flow rate at a given axial location. 

AIRMIX 

(3) 

Curve fit constants for dilution air entrainment. 

B 

(N) 

l+2y , 

BETA 

(N) 

Vector of intermediate coefficients £. 

EHRAT 


Burner diameter to height ratio. 

BNO 


Evaporation mass transfer number. 

C 

(N) 

-p-y , 

Cl 


Defined by equation 5-15. 

C2 


Defined by equation 5-16. 

C3 


Locally defined variable. 

CCONFG 

(2) 

Curve fit constants for conductivity of fuel vapor 
as f(T). 

CCPFG 

(2) 

Curve fit constants for specific heat of fuel vapor 
as f (T) . 

CHT 

<5) 

Curve fit constants for distance above orifice 
as f(TAU). 

CLAT 

(3) 

Curve fit constants for heat of vaporization as f(T) 

CONI 


Locally defined variable. 

C0N2 


Locally defined variable. 

CONG 


Conductivity of gas phase. 

CONL 


Conductivity of liquid fuel. 

CONST2 


Locally defined variable. 

C0NST3 


Locally defined variable. 

CPF 


Specific heat of liquid fuel. 

CPG 


Specific heat of gas phase. 

CVAP 

(5) 

Curve fit constants for vapor pressure as f(T). 

CVISFG 

<2> 

Curve fit constants for fuel vapor viscosity as f(T). 

D 

(N) 

Coefficient vector containing elements of TD. 

D1 


Locally defined variable. 

D2 


Locally defined variable. 

DELU 


U2-U1 . 

DIFF 


Diff usivity. 

DISAIR 


Dispersion air flow rate. 

DR 


Radial spacing. 

DT 


Ramped time increment. 

DTO 


Initial time increment. 

DTINC 


Time increment multiplier. 

EXTENT 


Extent of prevaporization. 

FI 

y 

Locally defined variable. 

F2 


Locally defined variable. 

FMW 


Molecular weight of fuel, . 

FNAME 


Chemical name of fuel. 

GAMMA 

(N) 

vector of intermediate coefficients y . 


FORTRAN Variable Definition 


HT 

HTVAP 

1 

IFP1 

JTIME 

L 

N 

N2MW 
NT I ME 

02 MW 
P 

PRNO 

PVAPOR 

R 

REDHT 

REYNO 

RHOG 

RHOL 

RS 

RSO 

RSP 

RU 

SYOW 

T 

TO 

TAU 

TBOIL 

TD 

TDROPO 

TDROP 

TEVAP 

THETAM 

TOTAIR 

TSURF 

U1 

U2 

V 

VISCG 

XF 

XN2 

X02 

XPLUS 

YF 

YFLG 

YFLG1 

YFS 

YN2 

Y02 


Burner height. 

Heat of vaporization. 

Local index. 

Locally defined variable. 

Time index. 

Local index. 

No. of radial meshes. 

Molecular weight of N2. 

No. of time steps. 

Molecular weight of 02. 

Total pressure. 

Prandtl No, 

Vapor pressure. 

Radial position. 

Height of the reducing sect/ or . 

Reynolds No. 

Gas phase density. 

Liquid fuel density. 

Droplet radius. 

Initial droplet radius. 

=RS . 

Universal gas constant. 

Inverse average molecular weight. 

Gas phase temperature. 

Initial gas phase temperature. 

Cumulative evaporation time. 

Boiling temperature of fuel. 

(N) Droplet temperature. 

Initial droplet temperature. 

Average droplet temperature. 

Total evaporation time. 

Non dimensional temperature evaluated by , 
interpolation. 

Total air flow rate. 

Droplet surface temperature. 

Gas phase velocity at time=TAU-DT. 

Gas phase velocity at time=TAU. 

(N) Vector containing the computed solution of TD. 

Gas phase viscosity. 

Fuel vapor mole fraction. 

N2 mole fraction. 

02 mole fraction. 

Locally defined variable. 

Fuel mass fraction. 

Liquid fuel mass fraction. 

Initial mass fraction of liquid fuel. 

Mass fraction of fuel at droplet surface. 

N2 mass fraction. 

02 mass fraction. 


PROGRAM VAPOR 


C 

C AAAAAAAAAAAAAAA^jAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

C AA SPRAY VAPORIZATION PROGRAM AA 

C *a THIS PROGRAM USES A DROPLET TRACKING TECHNIQUE TO AA 

C AA CALCULATE THE PREVAPORIZATION PARAMETERS ALONG THE AA 

C AA SPRAY PATH AA 

C AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

REAL N2MW 

COMMON FMW , N2MW , 02MW , CFG , CONG , DIFF , RHOG , VISCG 
COMMON RU,P,XF,X02,XN2,SY0W 
COMMON CCPFG( 2) ,CC0NFG(2) ,CVISFG(2) 

DIMENSION TD(51),A(51),B(51),C(51) ,D(51) 

DIMENSION CVAP ( 5 ) , CLAT ( 3 ) , AIRMIX( 3 ) , CHT( 5 ) 

CHARACTER AH FNAME 
C 

C INITIAL CONDITIONS — * 

C 

N=20 
NTIME=60 
TAU=0 . 0 

PRINT A,' ENTER DROPLET DIAMETER (MICRONS) 

READ (A, A) RSO 

RS0=RS0/2.0A1.0E-6 

RSP=RS0 

RS=RS0 

PRINT A , ' ENTER AIR TEMPERATURE (C) 

READ (A, A) TO 
T0=T0+273 . 15 

PRINT A,' ADIABATIC VAPORIZATION (1=Y ,0=N) 

READ <A r A) ADIAB 

PRINT A,' ENTER THE HEIGHT OF THE REDUCING SECTIONS (CM) : ' 
READ (A, A) REDHT 
READ (8,200) FNAME 

READ (8,210) FMW , RHOL , TBO IL , CONL , CPF , YFLG1 , TEVAP , BHRAT 
READ (8,220) (C’VAP( I) ,1 = 1,5) 

READ (8,230) (CLAT( I) ,1=1,3) 

READ (8,240) TOTAIR,DISAIR, (AIRMIX( I) ,1=1,3) 

READ (8,250) (CCPFG( I) ,1=1,2) 

READ (8,260) (CCONFG( I ) ,1=1,2) 

READ (8,270) (CVISFG(I) ,1=1,2) 

READ (8,280) (CHT( I) ,1=1,5) 

200 FORMAT (All) 

210 FORMAT (8F8.4) 

220 FORMAT (5E10.4) 

230 FORMAT (F8.1,F9.5,F10.5) 

240 FORMAT ( 2F7 . 1 ,F10. 1 ,2F7 . 1 ) 

250 FORMAT (2F9.6) 


260 FORMAT (2E10.3) 

270 FORMAT (2E10.4) 

280 FORMAT (5F8.4) 

02MW=32 . 0 
N2MW=28 . 0 
P=1.0 

RU*0. 08205 
T=TO 

TSURF=297. 15 
TDR0P=297 . 15 
TDR0P0=297. 15 
YF=0 . 0 
Y02=0. 231 
YN2=0 . 769 
DT0=2. OE-4 

DTINC= ( 2 . 0*TEVAP/NTIME-2. 0*DT0 ) / ( NTIME-1 ) 

Ul-27.59 
DO 5 1*1, N+l 
5 TD( I ) =TDROP 

YFLG=YFLG1*T0TAIR/DISAIR 
WRITE ( 5,101 > FNAME 

101 FORMAT ( 18X,A11 , / ) 

WRITE (5,102) 

102 FORMAT (' TIME(SEC) HT(CM) D(MIC) Tg(C) Ts(C) Td(C) 


* E' ) 

C 

C - BEGIN ITTERATION 

C 


DO 16 JTIME=1,NTIME 
DT=DTO+ ( JTIME-1 ) *DTINC 
TAU=TAU+DT 

IF (TAU.GE. 2 . 7E-4 ) GOTO 7 
AIRFLO=DISAJR+AIRMIX< 1 ) *TAU 
GO TO 9 


. 7 AIRFLO =AIRMIX ( 2 ) +AIRMIX ( 3 ) *TAU 

9 YFLG=YFLG*DISAIR /AIRFLO 

YF=YF*DISAIR/AIRFLO 

T= ( DISAIRAT+ ( AIRFLO-DISAIR ) *T0 ) /AIRFLO 
DISAIR* AIRFLO 
C 

C VELOCITY — 

C 


HT=CHT( 1 ) 

DO 18 1=2,5 

HT=HT+CHT(I)*(ALOG(TAU) ) ( 1-1 ) 

18 CONTINUE 

U2=( AIRFLO/ 6000 • 0) / ( 3 * 14159/4 . 0MHT/BHRAT)**2 ) 
HT=HT/ 100.0 


DELU=U2-U1 

C — VAPOR PRESSURE CALCULATIONS 

17 PVAPOR =CVAP ( 1 ) 

DO 20 1=2,5 

PVAPOR=PVAPOR+CVAP ( I ) A ( TSURF- 273.15)**(I-i> 

20 CONI' INUE 


non non o odd non 
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GAS PROPERTIES 


CALL FRACT <T,YF, YN2 , Y02) 

CALL PROPTY (T,YF,YN2,Y02) 

-- HEAT OF VAPORIZATION 

♦ 


HTVAP=CLAT( 1 ) 

DO 22 1=2,3 

HTVAP =HTVAP+CLAT ( I ) ATSURFAA ( 1-1 > 

22 CONTINUE 

YFS=1 . 0/ ( 1. 0- ( 1. 0-EXP (HTVAPAFMW/ 8315. 0A( 1 . 0/TSURF-1 . 0/TB0IL) ) > 

A / ( SYOWAFMW) ) 

IF (YFS.GT.YF) GOTO 40 
YFS=YF 

40 C0NST2=CLAT( 2 ) +2 . 0*CLAT< 3 ) ATSURF 

C0NST3 =EXP( HTVAPAFMW/ 3315 . 0 A ( 1 . 0 /TSURF-1 . 0 /TBOIL ) ) / ( SYOWAFMW ) A 
A (HTVAPAFMW/8315.0)/(TSURFATSURF)AYFSAYFS 

D2 = - ( CPGA ( T-TSURF ) A( YFS-1. 0) - < HTVAP+CPFA ( TSURF-TDROP ) ) A( YF-YFS) ) * 

A / ( CPGA < - ( YFS - 1 . 0 ) +CONST3 a ( T-TSURF ) ) +C0NST3 A ( HTVAP+CPFA 
A ( TSURF-TDROP ) ) - ( YF-YFS ) A ( CONST2+CPF ) ) 

TSURF=TSURF+D2 

IF (ABS(D2) .GE. 0.005) GOTO 17 
BNO = ( YF-YFS )/( YFS-1 . 0 ) 


DROPLET SIZE CALCULATIONS 

C1=0.25ASQRT(ABS(8.0A(9.81A(RHOG-RHOL)+DELU/DTARHOL)ARHOG • 

A /(18.0AVISCGAVISCG) ) ) 

C2=C0NG/ < CPGARHOL ) AAL0G< 1 . O+BNO ) 

C3 = 0 . 0 

DO 24 1=1,35 

C3=C3+(-1.0)AA(I+1)A(C1ARSAA1.5)AA(I-1)ARSARS/(1.5A(I-1)-»-2.0) f 

24 CONTINUE f 

38 FI =0.0 r ' ^ 

F2 = 0 . 0 
DO 30 1=1,35 

Fl=Fl4-( -1,0) A A ( 1+1 ) A ( C1ARSAA1 ,5)AA(I-1) ARSARS/ (1.5A(I-l)+2.0) 

F2=F2+( -1 . 0 ) aa < i+i ) A ( C1ARSAAI . 5 ) aa ( i-l ) ARS 
30 CONTINUE 

F1=C3-F1-C2ADT 

F2=-F2-C2 

D1=-F1/F2 

RS=RS+D1 

IF (ABS(Dl) .GE.1.0E-8) GOTO 38 
EXTENT= 1 . 0 - ( RS /RSP ) a A3 
U1=U2 

DROPLET TEMPERATURE I 

i 

TDROPO =TDR0P 

TD(N+1 ) =TSURF } 


i*) 
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DR=1 . 0/N 

A1=DTAC0NL/ (RHOLACFFARSARSADR) 

R=0 . 0 

DO 55 1=2, N 
R=R+DR 

A(I)=A1/R-A1AR/DR 
B( I ) =1 . 0+2. 0AA1AR/DR 
C( I ) =“A1/R-A1*R/DR 
D< I ) =TD( I ) 

55 CONTINUE 

C(2)=-2.0*A1 

D(N) =D(N) +(A1/ ( 1. 0-DR)+AlA( 1.0-DR) /DR) ATD(N+1) 

C 

C SOLVE A SYSTEM OF LINEAR EQNS TO GET DROPLET TEMPERATURE PROI 

C 

CALL TRIDAG ( 2 ,N,A,B,C,D,TD) 

TD( 1 ) =TD( 2 ) 

TDR0P=0 . 0 
DO 52 1 = 1, N 

C0N1= (TD( 14-1 ) -TD( I ) ) / 2 . 0 
C0N2=TD( 1+1 > -CONI A (FLOAT ( I ) ADR) 

TDR0P=TDR0P+3 . OA ( CONI /4 . OA ( ( FLOAT ( I ) *DR ) A*4- ( FLOAT ( 1-1 ) *DR) A> 
A +C0N2/3. 0A( <FL0ATU)ADR)A*3-(FL0AT(I-1)ADR)A*3) ) 

52 CONTINUE 
C 

C GAS PHASE TEMPERATURE 

C 

T= ( YFLGA ( CPF A ( TDROPO - ( 1 . 0- EXTENT ) ATDROP ) - EXTENT A ( HTVAP+TSURFA 
A ( CPF-CPG ) ) ) +CPGAT ) / ( CPGA ( 1 . O+EXTENTAYFLG ) ) 

IF (T.LE.TO) GOTO 14 
T=TO 

14 YF= ( YF+EXTENTAYFLG ) / ( 1 . 0 +EXTENTA YFLG ) 

YFLG= ( 1 . O-EXTENT) +YFLG/ ( 1 . 0 4- EXTENT A YFLG ) 

YN2=( 1 . 0-YF) / ( 1+0 .231/0.769) 

Y02=l . 0- ( YF+YN2 ) 

RSP=RS 

C 

C - THERMAL ENTRY LENGTH SOLN FOR WALL HEATING 

C 

IF ( ADIAB.EQ ,1.0) GOTO 11 

IF (HT.LT.REDHTA1.0E-2) GOTO 11 

PRNO=VISCGACPG/CONG 

REYNO =RHOG AU2 AHT7 BHRAT / V I SCG 

XPLUS = ( HT-REDHT* 1 . 0E- 2 ) / ( 1 . 43/200. 0 ) / ( REYNO APRNO ) 

IF (XPLUS. GT. 0.001) GO TO 62 
THETAM=1 . 0-XPLUSA ( 1 . 0-0 . 962 ) / 0 . 001 
GO TO 69’ 

62 IF (XPLUS. GT. 0.004) GO TO 63 

THETAM=0 . 962+( 0 . 001-XPLUS ) A ( 0 . 962-0.90$ ) / ( 0 . 004-0 . 001 ) 

GO TO 69 

63 IF ( XPLUS . GT .0.01) GO TO 64 

THETAM=0 . 908+ ( 0 . 004-XPLUS ) A ( 0 . 908-0 . 837 )/( 0 . 01-0 . 004 ) 

GO TO 69 

64 IF ( XPLUS . GT .0.04) GO TO 65 


THETAM=0.837+(0.01-XPLUS)A( 0.837-0. 628)/ (0.04-0. 01) 

GO TO 69 

65 IF (XPLUS.GT.0.08) GO TO 66 

THETAM^O .628+< 0 . 04-XPLUS ) * ( 0 .628-0 . 459 ) / ( 0 . 08-0 . 04 ) 

GO TO 69 

66 IF (XPLUS.GT.O.l) GO TO 67 

THETAM=0 . 459+ ( 0 . 08-XPLUS ) A ( 0 . 459-0 . 396 ) / ( 0 . 1-0 . 08 ) 

GO TO 69 

67 IF (XPLUS.GT.0.2) GO TO 68 
THETAM=0.369+(0.1-XPLUS)A< 0.396-0. 190) /(O. 2-0.1) 

GO TO 69 

68 THETAM=0. 0 

69 T=TO-(TO-T)*THETAM 
C 

11 WRITE (5,290) TAU,HT*100 . 0 , 2 . 0E6*RS,T-273 . 15 ,TSURF-273 . 15, 

A TDROP-273.15,BNO,1.0-(RS/RSO)A*3 
290 FORMAT <F8. 5,2X,F5. 2 ,F9 . 3 ,3 <F6. 2) ,F8. 5 ,F7 . 4 > 

16 CONTINUE 

STOP 
END 

SUBROUTINE FRACT(T,YF,YN2 ,Y02) 

REAL N2MW 

COMMON FMW , N2MW , 0 2MW , CPG , CONG , D IFF , RHOG , V I SCG 
COMMON RU,P,XF,X02,XN2,SY0W 
C 

C A A A A A A AAAAAA A AAAAA A A AAAAAA A A A A A A AAAA A A AAA A A A A A A A A A AAAA A A A A A A A>'; A 

C A MASS TO MOLE FRACTION CONVERSION AND AVERAGE DENSITY CALC. A 

C AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

SY0W=YF/FMW+Y02/02MW+YN2/N2MW 
XF=YF/ (FMWASYOW) 

X02 =Y02 i ( 02MWASY0W ) 

XN2 = YN2 / ( N2MWASY0W ) 

RHOG=P/(SYOWARUAT) 

RETURN 

END 


SUBROUTINE PR0PTY(T,YF,YN2,Y02) 

REAL N2MW 

COMMON FMW , N2MW , 02MW , CPG , CONG , DIFF , RHOG , VISCG 
COMMON RU,P,XF,X02 ,XN2 ,SYOW 
COMMON CCPFG ( 2 ) , CCONFG ( 2 ) , CVISFG ( 2 > 

C 


C A A A A A A A A A A A A A A A A AAAA AAAAA A A AAA A A AAA A A A A AAAAAAAAAAA AAA AAAAAA A A A A 

C + GAS PROPERTY CALCULATIONS , A 

C AAA AAA A A A AAAA A AAAAA A A A AAA A AAAAA AA AAAAAAAAAAAAAAAAAAAAAA A A AAAAAA 

C 

C CP MIXTURE 

C 


CPG= ( YFA ( CCPFG ( 1 ) +CCPFG ( 2 ) *T*1 . 8 ) +Y02* ( 0 . 36-5 . 375/SQRT(T*l . 8 ) 

A +47.8/(1.8AT) )+YN2A( 0.338-123. 8/(1. 8AT) +4. 14E4/(1.8A1.8ATAT) ) ) 
a A4186.9 


onoonoo n noon no 
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THERMAL CONDUCTIVITY OF MIXTURE 

CONG* ( XFA ( CCONFG ( 1 ) -*-CCONFG ( 2 ) AT4-CC0NFG ( 3 ) ATAT ) AFMWAA ( 1 . 0 / 3 . 0 ) 

A +X02A ( 0 . 00784914-6. 784E-5AT-4. 903E-9 ATAT) A02MWA A ( 1 • 0/ 3 . 0 ) 

A +XN2A ( 0 . 01439064-4 . 545E-5AT-H . 515E-9ATAT) AN2MWAA( 1 . 0/3 . 0 ) ) / 

A ( XFAFMWAA { 1 . 0 / 3 . 0 ) 4-X02A02MWAA ( 1 . 0 / 3 . 0 ) 4-XN2AN2MWAA ( 1 . 0 / 3 . 0 ) ) 


ABSOLUTE VISCOSITY OF MIXTURE 

VISCG* ( XFA ( CVISFG ( 1 ) 4-CVISFG ( 2 ) AT ) ASQRT(FMW) 
a 4-X02A(1.697E-54-5.51E-8AT-1.0E-10ATAT)ASQRT(02MW) 
a 4-XN2A ( 1 . 697E-5+5 • 51E-8AT-1 . 0E-10 ATAT) ASQRT( N2MW ) ) / 

A (XFASQRT(FMW)4-X02ASQRT(02MW)4-XN2ASQRT(N2MW) ) 

DIFF=C0NG/ (RHOGACPG) 

RETURN 

END 

SUBROUTINE TRIDAG( IF,L,A,B,C,D,V) 

DIMENSION A(l) ,8(1) ,C( 1 ) r D(l) ,V( 1 ) ,BETA( 51 ) f GAMMA( 51 ) 

AAAAAAAAAAAAAA A AAA A A A A AAA A A A AAAAA AA AAAAAAAAAAAAAA AAAAAAAAAA 
A LINEAR SIMULTANEOUS SYSTEM OF EQNS SOLVER A 

a REFs "APPLIED NUMERICAL METHODS , " CARNAHAN ET AL. A 
A JOHN WILEY & SONS , 1976. A 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

BETAt IF) =B( IF) 

GAMMA ( IF) =D( IF) /BETA( IF) 

IFP1 = IF4-1 
DO 1 I=IFP1 ,L 

BETA( I ) =B( I ) -A( I ) AC( 1-1 ) /BETA( 1-1 ) 

1 GAMMA(I)=(D(I)-A(I)AGAMMA(I-1) )/BETA(I) 

V(L) =GAMMA(L) 

LAST =L- IF 
DO 2 K=1 ,LAST 
I=L-K 

2 V(I)=GAMMA(I)-C(I)AV(I4-1)/BETA(I) 

RETURN 

END 
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E.3 Sorav Combustion 


FORTRAN Variable Definition 


AIRO 

AIRMAS 

B 

BRATE1 

BRATE2 

BRATE3 

BURND 

CF 

CO 

CPDT 

DIST 

DR 

DRO 

DRI (N> 

DRLST 

E 

ELLAST 

ELMAS 

FRATE1 

FRATE2 

FRATE3 

FUO 

FUMAS 

J 

NBRACH 

NDROP 

NODEN 

NO MW 

PHIO 

PI 

PORUT2 

Q 

RO 

REACT! 

REACT2 

REACT3 

RHO 

RHODIF 

RHOS 

RVDRDT 

SDRLST 

SOXF 

SYOWBK 

T (N) 

TAU 

TENLOG 


Total oxidizer mass around a burning droplet. 
Cumulative oxidizer mass. 

Preexponential constant for fuel oxidation reaction 
Backward reaction rate constant for: 02+N2 = 20+N2 

Backward reaction rate constant for: 0+N2 = NO+N. 

Backward reaction rate constant for: N+02 = NO+O. 

Burner height. 

Stoichiometric molar coefficient for fuel. 
Stoichiometric molar coefficient for oxydizer. , 

JCp dT. 

Center-to-center droplet spacing, 2 t 
= DRI. 

First radial spacing. 

Ramped radial spacing. 

Last radial spacing. 

Activation energy for fuel oxidation reaction. 
Locally defined variable. 

Locally defined variable. 

Forward reaction rate constant for: 02+N2 - 2Q+N2 

Forward reaction rate constant for: 0+N2 = NO+N. 

Forward reaction rate constant for: N+02 = NO+O. 

Total mass of fuel. 

Locally defined variable. 

Local index. 

Burning rate correction factor equation flag. 

Number of droplets in a planar array. 

Droplet number density. 

Molecular weight of NO. 

Overall equivalence ratio. 

3.14159 

Locally defined variable. 

Heat of reaction. 

Universal gas constant. 

Rate of 0 formation from: 02+N2 = 20+N2. 

Rate of N formation from: 0+N2 = NO+N. 

Rate of NO formation from: N+02 = N0+0. 

Gas phase density. 

P x 

Density at droplet surface. 

Locally defined variable. 

Locally defined variable. 

Stoichiometric oxygen to fuel mass ratio. 

Previous inverse average molecular weight. 

Flame temperature. 

Cumulative reaction time. 

Log base 10 of nondimensional particle spacing. 
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FORTRAN Variable Definition 


TGAS 


Initial reaction temperature. 

TGASO 


Inlet gas phase temperature. 

TLO 


Boiling temperature of fuel. 

TNEW 

(N) 

Updated flame temperature. 

TOTFUI , 


Sum of cumulative (reacted and unreacted) fuel ma.* 

TOTMAd 


Locally defined variable. 

TOTMSO 


Combined mass of fuel and oxidizer. 

TOTOX 


Sum of cumulative (reacted and unreacted) oxidize* 
mass . 

V 

(N) 

Radial velocity. 

XI 


Local x-coordinate of a droplet. 

XN 


Mole fraction of N. 

XNO 


Mole fraction of NO. 

XO 


Mole fraction of 0. 

Y1 


Local Y-coordinate of a droplet. 

YFGAS 


Initial fuel vapor mass fraction in gas phase. 

YFNEW 

(N) 

Updated fuel vapor mass fraction in gas phase. 

YN2 

(N) 

Mass fraction of N2 in gas phase. 

YN2GAS 


Initial mass fraction of N2 in the gas phase. 

YN2NEW 

(N) 

Updated N2 mass fraction in gas phase. 

YNO 

(N) 

Mass fraction of NO in gas phase. 

YNONEW 

(N> 

Updated NO mass fraction in gas phase. 

Y02 

(N) 

Mass fraction of 02 in gas phase. 

Y02GAS 


Initial mass fraction of 02 in the gas phase. 

Y02NEW 

(N) 

Updated 02 mass fraction in gas phase. 

YPMINS 


Combined mass fraction of product (no N2) 
at R+DR( 1-1 ) . 

YPPLUS 


Combined mass fraction of product (no N2) 
at R+DR ( I ) . 

YPROD 


Combined mass fraction of product (no N2) 
at R+DR(N) . 


141 


M 


PROGRAM COMB 
C 

C AAAA AAAAAAAAA A A AAAAAAAAAAAAAAAAAAAAA AAA A A A A A A A A A AAA A AAA A A A 

C * A SPRAY COMBUSTION PROGRAM A A 

C aa WITH NO FORMATION KINETICS AA 

C AAA AAAAAA A A A AAAAAAAAA AAA A A A AAAAAA A A A A A A A AAA A A A A A A A A A A A A A A A 

C 

REAL N2MW,N0MW 

COMMON FMW , N 2MW , 02MW , NOMW , CPG , CONG , DIFF , RHO , RU , P 
COMMON XF,X02,XN2,XN,X0,XN0,SY0W 

DIMENSION T( 301 ) ,TNEW( 301 ) , Yr’( 301 ) , YFNEW( 301 ) , Y02( 301 ) 
DIMENSION Y02NEWC 301 ) ,YN2( 301 ) ,DRI ( 301 ) 

DIMENSION YN2NEW<301) ,YN0(301) ,YN0NEW(301) ,V(301) 

C 

C ' INITIAL CONDITIONS - 

C 

351 FORMAT (F7.2,I5,F7.2,F5.2,F7.4> 

352 FORMAT (F7 . 2 ,F7 . 4 ,F6 . 1 ,F6 . 1 ,F6 . 2 ) 

353 FORMAT (4E10.3) 

READ ( 7 ,351 ) RS 0 , NODEN , TGAS 0 , PHI 0 , EXTENT 
READ (7,352) FMW,CONL,RHOL,TLO ,CO 
READ (7,353) HTVAP,Q,B,E 
TAU=0. 0 

RS0=RS0/2. OA1. OE-4 
TGAS0=TGAS0+273 . 15 
DR0-0.2 
CF=1 . 0 
02MW=32 . 0 
N2MW=28.0 
NOMW=30 . 0 
P=1.0 

RO =8. 315A 1000.0 
RU-0 . 08205 

S0XF=C0A02MW/ ( FMWAPHIO ) 

C 

C DROPLET CONFIGURATION IN THE ARRAY 

C 

CALL CONFIG (RSO,DIST, NODEN, NDROP) 

C 

RS0=RS0A1.0E-2 

RS=RSO 

N=INT( -DRO+O . 03+SQRT( (DRO-O . 03 ) AA2-4 . 0A( 1 . O-DIST) AO . 03 ) / 0 . 06 ) -1 
WRITE (5,35) DIST, 2 . 0 ARSO A 1E6 , NODEN, NDROP, PHI 0 
35 FORMAT (' INIT DROP SPAC. = ' ,F5. 2, ' DO «',F6.2,' NO. DENS 
A = ' , 14 , ' DROPS /ARRAY= ' ,13, ' PHI=',F5.3/) 

C 

FU0=4 . 0/3 . 0A3 . 14159ARHOLARSOAA3/ ( 1 . O-EXTENT) 

TOTMSO=FUOA ( 1 . O+SOXF/PHIO ) 

AIRO=FUOASOXF/PHIO 
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TGAS*I100. 0 

18 YFGAS =1 . 0/ (1 . 0+4 . 29ASOXF/ ( EXTENTAPHIO ) ) 

Y02GAS- < 1 « 0-YFGAS ) / ( 4 . 29 ) 

YN2GAS-3.29AY02GAS 

CPI/r=<YO2GASA{0.36A1.8A(TGAS-TGAS0>-2.0A5.375A(SQRT(TGASA1.8)- 
A 3QRT(TGAS0A1 . 8 ) ) +47 . 8A(AL0G( 1 . 8ATGAS ) -ALOG( 1 . 8ATGAS0 ) ) ) + 

A YN2GASA ( 0 . 338A1 . BA ( TGAS-TGASO ) - 123 . BA ( ALOG ( 1 . 8ATGAS ) - 
a ALOG( 1 . 8ATGAS0 ) ) -< 4 . 14E4 ) / 1 . 8A ( 1 , O/TGAS-1 . 0/TGASO ) ) +YFGASA 
A ( 0.0694A1.8A( TGAS-TGASO) 4- (5. 27E-4)/ 2. 0A1.8A1.8A( TGASATGAS- 
A TGASOATGASO) ) )A2326.0 
YFGAS ■ YFGAS -CPDT / Q 
Y02GAS =Y02GAS -CPDT/ QASOXF 
IF ( YFGAS. GT. 0.001) GOTO 23 
TGAS-TGAS-50.0 
GO TO 18 

23 DO 10 1=1, N+l 
T ( I ) -TGAS 
YF( I ) = YFGAS 
Y02( I ) =Y02GAS 
YN2( I ) =YN2GAS 
YNO( I ) =0 . 0 
V<I)*0.0 
10 CONTINUE 
T( 1 ) =TL0 
YF( 1 ) *0 . 80 
YC2(1)=0.0466 
YN2(1) =0.1534 
TNEW( 1 ) =T< 1 ) 

YFNEW(1)=YF(1) 

Y02NEW( 1 ) *Y02 ( 1 ) 

YN2NEW( 1 ) =YN2 ( 1 ) 

DT-RS0A0.7E-3 

DRLST-DR0+ ( FLOAT (N-l) )A0.06 
SDRLST=DRLST 
DO 20 1=2,301 

DRI ( I) =DRO+(FLOAT( 1-2 ) >ao.06 
20 CONTINUE 
C 

C BEGIN I ITERATION 

C 

DO 16 JTIME=1, 150000 
IF (JTIME.LT. 10) GO TO 43 
DT=800. OARSARS 
43 TAU=TAU+DT 

CALL FRACT (T,l ,YF,YN2 ,Y02 ,YNO) 

SYOWBK=SYOW 

RHOS=RHO 

C 

C — — DROPLET* INTERACTIONS — 

C 

CALL CONFAC(DIST,RSO f RS ,NDROP,ETA) 



C 

C 

C 


BEGIN RADIAL MARCHING 


nnnnn on non non n n n non 


IF (RS.GE. 1.5E-6) GOTO 339 

FUMAS=0.0 

GO TO 354 

339 FUMAS=4 . 0/ 3 . 0A3 . 14159ARH0LARSAA3 

354 TOTMAS *FUMAS 

AIRMAS=0. 0 
R=1 . 0 

DO 12 1=2, N 
DR=DRI( I ) 

R=R+DR 


GAS PROPERTIES — 

CALL FRACT <T,I ,YF,YN2 ,Y02 ,YN0) 

CALL PROPTY <T,I ,YF,YN2 , Y02 , YNO ) 

— — CONSERVATION OF ENERGY - 

P0RUT2=(P/<RUAT(I) ) >AA2 

TNEW( I ) =T ( I ) +DTA ( - V< I W ( RS ADR > A < T ( 14-1 > -T( I ~ 1 ) ) / 2 . 04-CONG/ ( RHOA 
A CPGARSARS)A( (T(I+1>-T(I-1) > / (RADR) 4-(T( 14-1 ) -2 . 0AT< I ) 4T( 1-1 > ) 

A / ( DRADR ) ) +CFAFMWAP0RUT2 / ( RHOACPG ) AXFAX02 ABAEXP ( -S/ ( RO AT( I) ) ) AQ) 

CONSERVATION OF FUEL - 

P0RUT2 = ( P / ( RU ATNEW ( I ) ) )AA2 

YFNEW< I ) =YF( I ) 4-DTA < —V ( I ) / ( RSADR ) A( YF( 14-1 ) -YF( 1-1 ) ) /2 . O+RHOADIFF/ 
A (RHOARSARS)A( (YF(I+1)-YF(I-1) )/(RADR)4-(YF(I+l)-2.0AYF(I)+ 

A YF( 1-1 ) ) / (DRADR) ) -CFAFMW/ RH0AP0RUT2 AXFAX02 A 
* BAEXP(-E/(ROATNEW(I) ) ) ) 

IF (YFNEW(I) .GE.0.0) GOTO 300 
YFNEW( I ) =0 . 0 


CONSERVATION OF OXIDIZER 

00 Y02NEW( I ) =Y02 ( I ) +DTA ( -V( I ) / ( RSADR ) A ( Y02 ( 1 4-1 ) -Y02 ( 1-1 ) ) / 2 . 04-RHOA 
A DIFF/(RHOARSARS)A( ( Y02( 14-1 ) -Y02 ( 1-1 ) ) / (RADR) +( Y02 ( 14-1 ) 

A -2.0AY02(I)+Y02(I-1) >/(DRADR) )-C0A02MW/RH0AP0RUT2AXFAX02A 
A BAEXP(-E/(ROATNEW<I)))> 

IF <Y02NEW( I) .GE.0.0) GOTO 400 
Y02NEW(I)=0.0 

444444444444444444444444444444444444444444444444444444444444444 

\ 

4 - NO KINETICS 4 

44444444444444444444444 + 4444444444444444444444444444444444 + 444+4 

- REACTION RATES FOR ZEFDOVICH MECHANISM 

400 FRATE1=6.0E10A4.94E8/(R0ATNEW(I))AEXP(“4.94E8/(R0ATNEW(I) ) ) 

BRATE1=1 . 0E8 

FRATE2=1.36E11AEXP(-3.15E8/(R0ATNEW(I) ) ) 
BRATE2=3.1E10AEXP(-1.4E6/(R0ATNEW(I) ) > 
FRATE3=6.4E6ATNEW(I)AEXP(-?,.62E7/(R0ATNEW(I) ) ) 
BRATE3=1.5E6ATNEW(I)ASXP< -1.62E8/ (RQATNEW I ) ) ) 


r W'f 1 "V* * 
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X0«SQRT(FRATE1/BRATE1*RUATNEW< I) /PAX02 ) 

XN=XOA( BRATE3AXNO+FRATE2AXN2 ) / ( FRATE3 AX02+BRATE2AXNQ ) 

C 

C SOURCE AND SINK TERMS FOR ZELDOVICH MECHANISM 

C 

REACT 1 “P0RUT2 A ( 2 . 0AFRATE1AXO2AXN2-2. 0ABRATE1AX0AX0AXN2A 
A SQRT( PORUT2 ) > 

REACT2 =PORUT2 A ( FRATE2 AXO AXN2 - BRATE2 *XNO *XN ) 

REACT3 *PORUT2 A ( FRATE3 AXN AXO 2 « BRATE3 AXNO AXO > 

C 

C - N2 CONSERVATION 

C 

YN2NEW< I ) *YN2 ( I ) +DT* ( -V ( I ) / ( RSADR ) A ( YN2 ( 1+1 ) -YN2 ( I - 1 ) ) / 2 . 0+RH( 
A DIFF/ ( RHOARSARS )*( ( Y,N2(I+1 ) -YN2 ( 1-1 ) ) / ( RADR ) + ( YN2 ( 1+1 ) 

A -2.0*YN2(I)+YN2(I-1) )/(DRADR) ) ) 

C 

C — - NO CONSERVATION — - — — 

C 

YNONEW< I ) «YNO ( I ) +DTA ( -V( I ) / ( RSADR ) A ( YNO ( 1+1 ) -YNO < 1-1 ) ) / 2 . 0+P.HC 
A DIFF/ ( RHOARSARS ) A ( ( YNO ( 1+1 ) -YNO (I-l))/( RADR ) + ( YNO ( 1+1 ) 

A - 2 . 0 AYNO < I ) +YNO (I - 1 ) > / ( DR ADR ) ) +NOMW/RHOA ( REACT2+REACT3 ) ) 

C 

C ----- - OVERALL O/F RATIO — 

C 

ELMAS*4.0/3*0A3.14159A<RAA3-<R-DR>AA3>ARSAA3AP/(RUA 
A ( SYOW+SYOWRK) /2 . OA(TNEW( I ) +TNEW( 1-1) ) / 2 . 0 ) 

YPMINS»1.0-(YFNEW(I-1)+Y02NEW(I-1>+YN2NEW<I-1) > 

YPPLUS-1 . 0-< YFNEW( I ) +Y02NEW( I ) +YN2NEW( I ) ) 

AIRMAS « ( ( YO 2NEW ( I > +Y02NEW < I - 1 ) ) / 2 . 0+ ( YPMINS + YPPLUS )/ 2 . 0 A 
A (1,0-114. /514 . 0) ) AELMAS+AIRMAS 

FUMAS*( <YFNEW(I)+YFNEW(I-1,) )/2.0+eYPMINS+¥PPLUS)/2.Ali4. /514. ) 
A AELMAS+FUMAS 

TOTMAS = ELMAS +TOTMAS 
SYOWBK=SYOW 
C 

C — - CONTINUITY — 

C 

IF (RS.LT.1.5E-6) GOTO 64 
RVDRDT s CONLA ( T ( 2 > -T ( 1 ) ) / ( RS ADROAHTVAP ) 

C V( I) =( (P/ ( RUATNEW< I ) ASYOW) -RHO ) / ( 3 . ADT) ARSA ( 1 . -RAA3 ) + 

C A RVDRDT)/(RHOARAR) 

V ( I ) =RVDRDTA ( RHOL-RHOS ) / ( RHO ARHOLARAR ) 

64 IF ( JTIME.NE. 1 . AND.FLOAT( JTIME) /1000. ,NE. JTIME/1000) GOTO 12 
IF (I.GT.2) GOTO 700 
WRITE (5,444) 

444 FORMAT ( ' R T V YF Y02 YN2 ' 

A , ' XN XO YNO') 

WRITE (5,123) R-DR,INT(T( 1 ) ) ,V( 1 ) ,YF( 1 ) , 

A Y02( 1 ) r YN2 ( 1 ) ,0.0,0. 0 ,YNO( 1 ) 

700 WRITE (5,123) R, INT(TNEW( I ) ) ,V( I ) ,YFNEW( I ) , 

A Y02NEW( I ) ,YN2NEW( I ) ,XN,XO,YNONEW( 1 ) 

123 FORMAT ( IX, F7 . 1 , IX, 15 ,F7 . 3 ,3F7 .4 ,3E10 . 3 ) 

12 CONTINUE 


, * 

R 
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------ ASSIGN THE NEW VALUES 


DO 14 I*2,N 
T ( I ) *TNEW ( I ) 

YF( I ) *YFNEW( I ) 

Y02 ( I ) =Y02NEW( I ) 

YN2< I ) ®YN2NEW( I ) 

YNO ( I ) =YNONEW( I ) 

14 CONTINUE 

IF (RS.LT.1.5E-6) GOTO 66 

o/F RATIO CHECK — — — 

ELLAST=TOTMSO-TOTMAS 

YPR0D= 1 . 0 - ( Y02NEW < N ) +YFNEW ( N ) +YN2NEW (N) ) 

72 TOTOX= AI RMAS +ELLAST A ( Y 0 2NEW ( N ) + YPRODA (1.0-114. /514. ) ) 

TOTFUL=FUMAS+ELLASTA ( YF(N) +YPR0DA114 . / 514 . ) 

IF <ABS(TOTOX/TOTFUL-SOXF) .LT.0.01) GOTO 74 
ELLAST= ( SOXFAFUMAS-AIRMAS ) / < Y02 ( N > +YPRODA (1.-114/514. 

A -S0XFA114./514. )-SOXFAYF(N) ) 

GO TO 72 

7 4 DRLST = ( ( ELLAST ) / 

A (4 . / 3 . A3 . 14159ARSAA3AP/ ( RUASYOWATNEW( N ) ) > +RAA3 ) AA ( l . 0/ 3 . 0 ) -R 
SDRLST =DRLST+SDRLST 
IF < SDRLST. GE. DRI (N+l ) ) GOTO 45 
N=N-1 

SDRLST =SDRLST-»-DR I ( N+l ) 

45 DRLST=SDRLST 

IF < SDRLST. LE. DRI ( N+2 ) ) GOTO 41 
SDRLST=SDRLST-DRI ( N+l ) 

N=N+1 

DRLST=SDRLST 

BOUNDARY CONDITIONS AT SURFACE 

41 RVDRDT=C0NLA(T(2)-T(1) )/(RSADROAHTVAP) 

RS =RS - DTARVDRDT / RHOL AETA 
IF <RS) 99,99,24 

24 V ( 1 ) -RVDRDTA (1.0/ RH0S+ , O/RHOL) 

RHODIF=CONL/ ( ( 0 . 0694+5 . 27E-4ATL0A1 . 8 ) A4186 . 9 ) 

YF ( 1 ) = ( RHODIFAYF ( 2 ) / ( RSADRO ) +RVDRDT ) / ( RHODIF/ ( RS ADRO ) +RVDRDT ) 
Y02 ( 1 ) = ( RH0DIFAY02 ( 2 ) / ( RSADRO ) ) / ( RHODIF/ ( RSADRO ) +RVDRDT ) 

YN2 ( 1 ) = ( RHODIFAYN2 ( 2 ) / ( RSADRO ) ) / ( RHODIF/ ( RSADRO ) +RVDRDT ) 

YNO ( 1 ) = ( RHODIFAYNO ( 2 ) / ( RSADRO ) ) / ( RHODIF/ ( RSADRO ) +RVDRDT ) 

GO TO 68 

66 DO 62 1 = 1 Ji 
V(I)=0.0' 

62 CONTINUE 

YF( 1 ) =YF( 2 { 

YN2 ( 1 ) =YN2 ( 2 ) 

Y02 ( 1 ) =Y02 ( 2 ) 

YNO( 1 ) =YNO( 2 ) 

C T(U=T(2) 
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BOUNDARY CONDITIONS AT AMBIENCE 


68 T(N+1)*T(N) 

7F(N+1)=YF(N) 

Y02(N+1) =Y02(N) 

YN2( N+l ) =YN2(N) 

YNO ( N+l ) =YN0 ( N ) 

V(N+1)=V(N) 

28 IF < JTIME.NE.l. AND. FLOAT* JTIME)/ 1000. .NE.JTIME/ 1000) GOTO 16 
WRITE (5,222) TAU,2.0ARSA1.OE6,1.O-(RS/RSO) AA3,ETA, 

A JTIME,DR,N 

222 FORMAT (' T(SEC) *' ,E10.3, ' ;D(MIC)=' ,F6.3, ' ;E=' ,F5.3, 

A ' ;ETA= ' ,F5. 3 , ' ; ITER= ' , 16 , ' ?DEL= ' ,F5 . 1 , ' ;N= 1 ,13/) 

16 CONTINUE 
99 STOP 
END 

SUBROUTINE FRACT(T,I ,YF,YN2 ,Y02 ,YN0) 

REAL N2MW,N0MW 

COMMON FMW , N2MW , 0 2MW , N0MW , CPG , CONG , DIFF , RH0 , RU , P 
COMMON XF,X02,XN2,XN,X0,XN0,SY0W 
DIMENSION T ( 1 ) , YF ( 1 ) „YN2( 1 ) , Y02 Cl ) ,YNO( 1 ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAj 1 

A MASS TO MOLE FRACTION CONVERSION AND AVERAGE DENSITY CALC. 
AAA A A A A A A AAA AAA A A A AAA A A AAA AAA A A A AAA A A AAA AAA A A A AAA A A AAA AAA A A hi 

SY0W=YF( I ) /FMW+Y02 ( I ) /02MW+YN2 ( I ) /N2MW+ ( 1.0~YF( I) -Y02 < I ) -YN2 f 
A A( <114./ 514. )/ 18. 0+(l. -114. 0/514.0/44.0) 

XF=YF( I ) / (FMWASYOW) 

X02=Y02 (I ) / ( 02MWASY0W) 

XN2 = YN2 ( I ) / ( N2MWAS YOW ) 

XN0=YN0( I ) / (NOMWASYOW) 

RHO=P/(RUAT(I)ASYOW) 

RETURN 
END ' 

SUBROUTINE PROPTY ( T , I , YF , YN2 , Y02 , YNO ) 

REAL N2MW,N0MW 

COMMON FMW , N2MW , 02MW , NOMW , CPG , CONG , DIFF , RHO , RU , P 
COMMON XF,X02,XN2,XN,X0,XN0,SY0W 
DIMENSION T( 1 ) ,YF(1) ,YN2( 1) ,Y02( 1) ,YN0(1) 

A A A A AAA AAA A A A A A AAA AAA A A A A A AAA AAA A AAA AAA A AAA A AAA AAA A A A AAA AAA A A 

A GAS PROPERTY CALCULATIONS 

A A A A A A A A A A A A A A A A A A AAA AAA A A A A A A A A A A A A A A A A A AAA AAAA AAAA A A A A A A A A A 

CP MIXTURE — ■! 

CPG=(YF(I)A(0.0694+5.27E-4AT(I) *1 . 8 ) +Y02 < I ) A < 0 . 36-5 . 375/ 

* SQRT(T(I)*1.8)+47.8/Q.8*T(I) ) ) +YN2( I > A( 0 . 338-123 . 8/ ( 1 . 8*T< II 
A +4.14E4/(1.8*1.8*T(I)AT(I) ) ) +( 1 . 0-YF< I ) -Y02 ( I ) -YN2 ( I ) )A 
A ( 0. 681A ( 0.368-148. 4/ (1.8AT( I) ) +3 . 2E4/ < 1 . 8A1 . 8AT< I ) AT< I ) ) ) 


A +0.319M1.102-33.1/SQRT(1.8*T(I) )+416.0/(1.8*T(I) ) ) ) )*4186.9 

THERMAL CONDUCTIVITY OF MIXTURE — 

CONG=(XF* ( ~0 . 0268344-1 . 329E-4AT* I ) -1 . 741E-8AT< I ) at ( I ) ) A 
A FMWAA(1. 0/3. 0)+X02A(0. 0078491+6. 784E-5AT(I)-4.903E-9ATa)A 
A T( I ) ) A02MWAA ( 1 . 0/ 3 . 0 ) +XN2A ( 0 . 0143906+4 . 545E-5*T( I ) +1 . 515E-9A 
A T( I ) AT ( I ) ) AN2MWA A ( 1 . 0/3 . 0 ) + ( 1 . 0-XF-X02 -XN2 ) A ( 0 . 681 A 
A (-0.0009023+6.915E-5AT(I)-3.688E-9AT(I)AT(I) )A44.0 AA<i.o/ 3.0) 
A +0. 319A(0. 0048332+1. 303E-4AT(I)-7. 059E-9AT(I)AT(I) >A 
A 18. 0AA(l. 0/3.0) ) )/(XFAFMWAA(l.o/3.0)+X02A02MWAA(l.o/3.0)+XN2A 
A N2MWAA( 1.0/3.0)+(l. 0-XF-XN2-X02 ) 

A A(0.681A44.0AA(1.0/3.0)+0.319A18.0AA(1.0/3.0) ) ) 

DIFF=C0NG/ (RHOACPG) 

RETURN 

END 

SUBROUTINE CONFIG ( RSO ,DIST,NODEN,NDROP ) 

AAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAA AAA A AAA AAAAAAAA 

A PROGRAM TO DETERMINE DROPLET SPACINGS FOR AN EVENLY A 
a DISPERSED , MONOSIZED ARRAY A 

A AAA AAA AAA AA AA AAA A A AAA AAA AAA AAA A A AAA A A A AAA A A A A A A A A A AAA AAA A A A A 

NDROP=0 
BURND=1 . 43 
PI=3. 1415927 

DIST=(2.0/ (N0DEN*SIN(PI/3 . 0 ) ) )**( 1 . 0/3 . 0 ) / ( 2 . OARSO ) 
BURND=BURND/ ( RSO +2 . 0 ) 

DO 15 J=l,INT(BURND/(2.0ADISTASIN(PI/3.0> ) )+l 
DO 20 I=1,INT(BURND/(2.0ADIST) )+l 
X1=(I-1)ADIST 

IF ( FLOAT (J)/2.0.NE.INT(J/2)) GOTO 30 

Xl=Xl+DIST/2 . 0 

Y1 = ( J-l ) ADIST+SIN ( PI / 3 . 0 ) 

IF (SQRT(X1*X1+Y1*Y1) .GT.BURND/2 . 0) GOTO 20 
NDROP=NDROP+l 
IF (IAJ.EQ.l) GOTO 20 
IF(Yl.EQ.O.O) GOTO 40 
NDROP^NDROP+l 
IF (X1.EQ.0.0) GOTO 20 
NDROP=NDROP+2 
GO TO 20 
NDROP=NDROP+l 
CONTINUE 
CONTINUE 
RETURN 
END 
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SUBROUTINE CONFAC ( DIST , RSO , RS , NDROP , ETA ) 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
A SUBROUTINE TO CALCULATE BURNING RATE CORRECTION FACTOR A 
A FOR MONODISPERSE SPRAYS. STATEMENT NUMBERS CORRESPOND A 
A TO THE HT)f42ER OF INTERACTING DROPLETS IN A PLANAR ARRAY. A 
A A AAAAAA AA.-fcAAv- AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


TENLOG=ALOG10(2.0*DIST*RSO/RS) 

IF (TENLOG.LT* 3. 17) GOTO 25 
ETA =1.0 
GO TO 27 

25 NBRACH= (NDROP- 1 ) /6-8 

GO TO (55,61,67,73,79,85,91,97,103,109,115,121,127,133,139, 

A 145,151,157) ,NBRACH 

55 ETA=-1.2436435+TENL0G*< 2. 2164602-0. 8041338*TENLOG+0.1235272A 

A TENLOGATENLOG-O. 0064171ATENLOGAA3) 

GO TO 27 

61 ETA=-1.218263+TENL0G*( 2, 1265647-0. 7282704*TENL0G+0. 098952* 

A TENLOGATENLOG-O. 0036453 ATENLOGAA3) 

GO TO 27 

67 ETA=-1 . 2337979+TENLOG* (2.0958116-0. 6898761*TENLOG+0 . 0856408 A 

A TENL0GATENL0G-0.0021527ATENL0GAA3) 

GO TO 27 

73 ETA=-1.1997394+TENL0G*( 1. 9788084-0. 5967735*TENLOG+0. 0559048* 

A TENLOGATENLOG+O . 0012112ATENLOGAA3 ) 

GO TO 27 

79 ETA=-1 . 1988726+TENL0G*< 1 . 9448054-0. 5593372*TENLOG+0 . 0432551A 

A TENLOGATENLOG+O. 0026203ATENL0GA A3) 

GO TO 27 

85 ETA=-1.1423215+TENL0G*<1. 8058571-0. 4514555ATENLOG+0 . 0088748* 

A TENLOGATENLOG+O . 0065384ATENLOGAA3 ) 

GO TO 27 

91 ETA=-1 . 1428598+TENL0GA ( 1 . 7650879-0 . 411 5640 ATENLOG-O .0042657A 

a TENLOGATENLOG+O . 0079962ATENLOGAA3 ) 

GO TO 27 

97 ETA=-1 . 0798480+TENLOGA ( 1 . 6162612-0 . 2968254ATENLOG-0 . 0409864A 

a TENLOGATENLOG+O . 0122215ATENL0GAA3 ) 

GO TO 27 

103 ETA=-1 . 0748937+TENLOGA ( 1 . 5707712-0 . 2552938ATENLOG-0 . 0544528A 

A TENLOG ATENLOG+ 0.0137137 ATENLOG A A 3 ) 

GO TO 27 

109 ETA=-1 . 0080442+TENLOG* ( 1 . 4162596-0 . 1363881+TENLOG-O . 0927271A 

a TENLOGATENLOG+O. 01 81601 ATENLOGA A3) 

GO TO 27 

115 ETA=-0.9990754+TENLOG*(1. 3675382-0. 0938431*TENLOG-0. 1063837* 

A TENLOGATENLOG+O. 0196740ATENL0GAA3) 

GO TO 27 

121 ETA 25 -0 . 9380671+TENL0G* ( 1 . 2264606+0 . 0151967+TENL0G-0 • 1417152* 

A TENLOGATENLOG+O. 02 38 125 ATENLOG* *3) 

GO TO 27 
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127 ETA=-0 . 9125118+TENLOGA ( 1 . 14763254-0 . 0791626 ATENLOG-C .1622019* 

* TENLOGATENL0G+0.0261359ATENLOGAA3) 

GO TO 27 

133 ETA=-0. 850963 1+TENLOGA ( 1 . 006496+0. 1884841 ATENLOG- 0. 1978621 A 

A TENLOG ATENLOG+ 0.0303513 ATENLOG A A 3 ) 

GO TO 27 

139 ETA--0 . 8368802+TENL0GA( 0 . 9543458+0 . 23 17049 ATENLOG- 0 . 21 15844 A 

A TENL0GATENL0G+ 0.0318774 ATENLOG A A 3 ) 

GO TO 27 

145 ETA=-0.7636684+TENLOGA<0. 7898749+0. 3588799ATENL0G-0 . 2532952A 

A TENLOGATENLOG+O . 0368515ATENL0GAA3 ) 

GO TO 27 

151 ETA=-0 . 7478763+TENL0GA (0.7367701+0. 4021712ATENL0G-0 . 26 7 002 5 A 

A TENLOGATENLOG+O. 0383796ATENLOGAA3 ) 

GO TO 27 

157 ETA=-0. 6827373+TENLOGA (0.5896260+0. 5165399 ATENLOG-O . 3047329 A 

A TENLOGATENLOG+O. 04 2 9 093 ATENLOG A A3) 

27 RETURN 
END 
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E.4 Sample Input Data 

E. 4 . 1 Droplet Interactions Program 


/ 

6 


• 


0.000 

0.000 

0.000 

1.0 

20.000 

0.000 

0.000 

1.0 

-20.000 

0.000 

0.000 

1.0 

10.000 

17.321 

0.000 

1.0 

10.000 

-17,321 

0.000 

1.0 

-10.000 

-17.321 

0.000 

1.0 

-10.000 

17.321 

0.000 

1.0 


E.4. 2 Sprav Vaporization Procrram 


N- OCTANE 


114,2 707.0 398.15 

0.132 

2205.0 

0.0663 0.1703 12.45 

2.267E0 0.2138E0 8 

368370.0 421.2021 -1. 

. 594E-3 
4705 

1.42E-4 

1.28E-5 

3400.0 1300.0 7253000. 
0.0694 0.000527 

-1.306E-1 7.997E-4 

7 . 2029E-7 1 . 7494E-8 

0 3272.0 

498.0 


31.59 9.935 1.372 

0.095 

0.0026 



E . 4 . 3 Sprav combustion P ro crram 

47.68 400 23.64 1.00 0.43 

114.20 0.0199 707.0 380.0 12.50 

3.027E5 4.460E7 1.150E11 1.260E8 
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